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ABSTRACT 
Consider  the  large  systems  of  linear  equations  AhUh  =  fh  which  arise  from 
the  discretization  of  a  second-order  elliptic  boundary  value  problem.  Con- 
sider also  the  preconditioned  systems  (i)  S^"U/,u/,  =  Bh^fn;  and  (ii) 
AhBh^Vh  =  fh,  Uh  =  Bh^h  where  Bh  is  itself  a  matrix  which  arises  from  the 
discretization  of  another  elliptic  operator.  We  discuss  the  effect  of  boundary 
conditions  (of  A  and  B)  on  the  L2  and  Hi  condition  of  Bh^Ah,  AhBh^.  In 
particular,  in  the  case  of  H2  regularity  one  finds  that  |5a"U;,Il,  is  uniformly 
bounded  if  and  only  if  A*  and  B*  have  the  same  boundary  conditions  while 
IAaSa~Hl  is  uniformly  bounded  if  and  only  if  A  and  B  have  the  same  boun- 
dary conditions.  Similarly,  IS^'U^I/z^  is  uniformly  bounded  only  if  A  and  B 
have  homogeneous  Dirichlet  boundary  conditions  on  the  same  portion  of  the 
boundary.  This  latter  result  does  not  depend  upon  H2  regularity. 
Key  words,  preconditioning,  elliptic  operators,  condition  number,  Hi  regu- 
larity, boundary  conditions 

1.   Introduction 

Let  fl  C  R2  be  a  bounded  domain  with  a  boundary  dH.    Consider  a  uniformly  elliptic 
operator  A  of  the  form 

Au  =   -KanUx)x  +  (ai2Uy)x  +  (ai2Ux)y  +  (<222«y)y] 

(1.1a) 
■¥  aiUx  +  a2Uy  +  aQU 

with  boundary  conditions  of  the  form 

«  =  0   .     on     ro=  ro(A)  (Lib) 

1^  =  ao(tT)u(a)  +  ai((T)  -|^    ,     on     Fi  =  Y^iA)  (Lie) 

where  Yq{A),Y\{^)  are  a  partition  of  dYl,  i.e., 

an  =  Yq{A)  u  ri(A) 

and  -|^  denotes  the  exterior  normal  derivative  while  -r^  denotes  the  tangential  derivative 

dV  OCT 

and  a  denotes  the  tangential  variable.   The  uniform  ellipticity  of  A  is  expressed  by  two  posi- 
tive constants  0  <  X.(A)  ^  A(A)  such  that,  for  all  {x,y)  €  H  we  have  the  inequalities 

MA)  [42  +  Ti^]  <  aii42  +  2ai2§Ti  +  ^22^^  ^  A(A)[^2  +  ^2j    _  (12) 

Consider  a  discretization  A/,  of  this  elliptic  operator  and  the  system  of  linear  equations 

AhlL  =  l  (1.3a) 


which  arises  in  the  numerical  solution  of  the  boundary-value  problem 

Au  =  f   .  (1.3b) 

Our  concern  is  with  the  effective  solution  of  (1.3a).  In  certain  cases,  e.g.,  the  self-adjoint 
Laplace  operator  with  Dirichlet  boundary  conditions  (Fq  =  afl.ri  =  0)  in  a  square,  there 
are  relatively  fast  and  stable  methods  for  the  solution  of  this  problem.  However,  in  the  gen- 
eral case  of  nonself-adjoint  problems  whose  symmetric  part  may  be  indefinite,  this  can  be  an 
expensive  and  time  consuming  problem.  Part  of  the  difficulty  arises  because  the  condition 
number  of  the  discrete  approximations, 

ChiAn)  :  =  UJa  •  \Ah\    ,  (1.4) 

is  of  order  h~^.  Hence,  one  approach  has  been  to  "precondition"  Ah  to  obtain  a  problem 
with  a  smaller  condition  number  (cf  [BP],  [CoGu],  [DoDu],  [Dy],  [ElSh],  [Gu],  and  [Wi]). 
That  is,  given  an  appropriate  Bh  one  considers  either 

Bh^AhU.  =  Bh^l  (1.5a) 


or 

'.-1„  =  f        n.-l 


AhBh'^  =  i  ,  Bh^y.  =  u.   .  (1.5b) 

In  particular,  it  has  been  suggested  that  Bh  be  the  discretization  of  another  elliptic 
operator  B  which- has  been  chosen  so  that,  among  other  things,  Bh^  is  relatively  easy  to 
obtain.  In  addition,  it  is  frequently  suggested  that  B  be  chosen  as  the  "leading  part"  of  A 
with  the.  same,  boundary  conditions  aa  A..  However,  as  we  shall  see,  if  one  preconditions  on 
the  right,  i.e.,  (1.5b),  then  this  is  a  reasonable  recipe  while  if  one  preconditions  on  the  left, 
i.e.,  (1.5a)  then  it  is  better  to  choose  B  so  that  A*  and  B*  have  the  same  boundary  condi- 
tions. The  situation  is  further  complicated  by  the  question:  which  norm  should  be  used  in 
the  computation  of  the  condition  number?  While  a  careful  look  at  the  theory  of  iterative 
processes  suggests  that  the  L2  (or  I2)  norm  is  the  correct  choice,  much  of  the  literature  (for 
finite-element  equations)  deals  with  the  Hi  norm. 

The  goal  of  this  paper  is  to  clarify  this  situation.  For  this  purpose,  we  will  discuss  both 
the  L2  condition  numbers  Q/S^'U/,),  Ci^iAhBh^)  and  the  Hi  condition  numbers 
Cn.iBh^Ah).  However,  in  order  to  do  this  we  must  discuss  the  elliptic  partial  differential 
operator  and  the  discrete  approximation  with  greater  care  and  precision. 

Let  /  ^  L2(fi)  and  consider  the  boundary-value  problem 

Au  =  f   .  (1.6) 

Suppose  that,  in  some  sense,  this  problem  has  a  unique  solution  for  all  /  ^  L2(il).  That  is, 
we  suppose  that  in  some  general  sense  A  is  invertible  on  Z,2(ft).  If  the  coefficients  of  (1.1a), 
(1.1c)  are  smooth,  dD,  is  smooth  and  either  Tq  =  0  or  Ti  =  0  then  (see  [Or])  A~^  is  a 
bounded  map  from  LjiO.)  to  Hii^).    In  other  words,  there  is  a  constant  KiiA)  >  0  such 
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that,  for  every  given  f  ^  L2  the  solution  u  of  (1.6)  is  in  H2(^),  satisfies  the  boundary  cxindi- 
tions  (1.1b),  (1.1c)  in  the  trace  operator  sense  and 

\u\H^^KiiA)\f\L^   .  (L7) 

Assuming  these  conditions  are  met  we  may  describe  the  operator  A  as  follows:  the  domain 
of  A  is  given  by 

D(A)  =  {m  ^  Hii^):   u  satisfies  (1.1b),  (1.1c)}  (1.8) 

Of  course,  A  is  a  bounded  map  from  its  domain  onto  L2.  That  is,  there  is  a  constant  K2(A) 
such  that 

UmL,  ^  K2(A)\u\h^   ,    u  €  DiA)    .  (1.9) 

In  the  general  case  when  neither  Tq  =  0  or  Fi  =  0  the  situation  is  somewhat  delicate. 
For  example,  let  A  be  the  Laplace  operator.  Then  one  can  exhibit  (polygonal)  domains  and 
boundary  conditions  for  which  (1.7)  is  false-(see  [Gr]).  In  this  latter  case,  some  solutions 
(with  Au  =  /  €  L2)  are  not  in  //2(n).  Our  discussion  of  the  L2  condition  of  AhBh^  and 
Bh^Ah  depends  on  the  estimates  (1.7),  (1.9)  and  the  constants  Ki(A),K2iA),  KiiA*).  K2{A'), 
K\i,B),  K2{B),  KiiB'),  K2(B*).  Hence,  while  our  discussion  is  sometimes  phrased  in  the 
generality  of  boundary  conditions  (1.1b),  (1.1c),  the  reader  is  reminded  that  the  H2  regular- 
ity is  usually  essential  to  the  discussion.  Of  course,  this  includes  a  large  class  of  problems. 
On  the  other  hand,  our  discussion  of  the  Hi  condition  does  not  require  H2  regularity  and 
hence  includes  the  case  of  the  general  boundary  conditions  (1.1b),  (1.1c). 

In  discussing  the  discretizations.  A/,,  it  is  convenient  to  adopt  the  language  of  finite- 
element  theory.  However,  our  discussion  is  fairly  general.  Imagine  a  family  of  finite 
dimensional  function  spaces  {Sh}  indexed  by  a  discretization  parameter  /i  -  0.  Every  func- 
tion V;,  €  Sh  is  described  by  a  finite  vector  £  =  {vj  :  j  =  1,2,...,/!^}  depending  on  the  choice 
of  basis  for  5^.  The  operator  Af,  of  (1.3a)  acts  on  this  representing  vector  £.  The  norm  used 
in  (1.4)  is  given  by 

<li,yL>h  =  h'^l.UjVj  (1.10a) 

\Li\l  =  <lL,U>h    ,  (1.10b) 

where  d  =  dimension  H  (c?  =  2  in  our  current  discussion).  We  will  assume  that  these  norms 
are  uniformly  equivalent  to  the  L2  norms  of  the  functions  v.  That  is,  there  is  a  constant 
c  >  0  such  that 

For  this  reason,  while  we  are  concerned  with  ChiBh^Af,)  or  ChiAhBh^),  it  suffices  to  con- 
sider the  L2  condition  numbers  ClJ.Bh^Ah)  or  Ci^iAhBh^).    Of  course,  the  L2  norm  is  given 
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by 

if, 8)  =  IJfgdxdy    ,     \f\l^=  (fj)    .  (1,10c) 

Within  the  L2  theory  our  basic  results  are  easily  stated. 

I).     If  we  seek  an  invertible  elliptic  operator  B  so  that  we  obtain  an  estimate  of  the 
form 

CiJiAhBhh  =  \AHBir\,  ■  \BhAh\  ^K    ,     SiWh^ho   ;  (1.11a) 

then  it  is  necessary  (and  essentially  sufficient)  that 

DiA)  =  D(5)    .  (1.11b) 

II).    If  we  seek  an  invertible  elliptic  operator  B  so  that  we  obtain  an  estimate  of  the 
form 

CiJiBh^AH)  =  \Bh^AH\L,  ■  Ua-^sJl,  ^K    ,     aUh^ho   ■  (I.i2a) 

then  it  is  necessary  (and  essentially  sufficient)  that 

D(A')  =  D(fi')    .  (1.12b) 

Much  of  the  previous  analysis  of  these  preconditioning  strategies  has  been  based  on 
problems  with  Dirichlet  boundary  condition  on  all  of  aft,  i.e., 

To  =  aft   .     Fi  =  0   . 

Unfortunately,  from  the  point-of-view  of  insight  into  the  general  case,  the  Dirichlet  problem 
is  misleading  because,  in  that  case 

DiA)  =  D(A*)    . 

There  is  one  notable  exception,  namely,  the  work  of  Bramble  and  Pasciak  [BP]  who  proved 
a  basic  result  about  such  preconditionings  (see  Theorem  5.1).  Unfortunately  they  failed  to 
appreciate  (1.12a),  (1.12b).  Hence,  in  their  example  they  used  a  preconditioner  which  satis- 
fied (1.11b)  not  (1.12b).   Thus,  their  one  L2  example  is  incorrect. 

In  fact,  with  the  already  noted  exception  of  [BP],  most  of  the  earlier  work  on  such 
elliptic  preconditionings  has  not  only  been  based  on  the  Dirichlet  problem  but  has  been 
based  on  the  Hi  condition  number,  aoL  the  L2  condition  number.  Most  of  the  analysis  of  the 
Hi  condition  number  is  based  on  the  concept  of  "Spectral  Equivalence"  introduced  by 
D'Yakonov  [DY].  Let  {Ah},  {Bh}  be  families  of  positive  definite  self-adjoint  operators 
defined  in  {5/,}.  These  families  are  said  to  be  uniformly  "spectrally  equivalent"  if  there 
exist  constants  ci,C2,  independent  of  h  such  that 
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While  D'Yakonov  considered  only  the  finite-dimensional  case,  it  is  easy  to  extend  this  con- 
cept to  the  infinite  dimensional  case  (see  [FMP]).  Let  a(u,v),  biu,v)  be  two  positive  defin- 
ite bilinear  forms  defined  on  a  Hilbert  space  V.  We  say  that  a  and  b  are  spectrally 
equivalent  if 

^        b(u,u)  '■ 

When  dealing  with  the  /a  condition  number  the  appropriate  concept  is  Li  norm 
equivalence.  This  concept  was  (implicitly)  used  in  [BP].  A  systematic  study  of  norm 
equivalence  is  given  in  [FMP].  For  our  purposes  we  use  the  following  definitions.  Let  A 
and  B  be  two  invertible  elliptic  operators.  We  say  that  A  and  B  are  L2  norm  equivalent  on 
L2(ft)  if  there  is  a  dense  set  D  C  L^iX^  and  a  constant  c  >  0  such  that:  for  every  f  i^  0, 
f  ^  D  we  have 

Ub-I/Iz.,  ^  c\f\L^,  |5A-1A^  <  c\f\L^   .  (1.13) 

We  say  that  A"^  and  B~^  are  L2  norm  equivalent  on  LiiO.)  if  there  is  a  dense  set  D  C  L2(fi) 
and  a  constant  c  >  0  such  that:   for  every  u  ^  D  we  have 

]A-^Bu\l^<  c\u\l^   ,     |5-U«L^  <  clwltj  (1.14) 

Similarly,  let  {A^}  and  {Bf,}  be  families  of  operators  defined  on  {Sh}  (presumably  Ah  and  Bh 
are  discretizations  of  A  and  5  respectively).  We  say  that  {Ah},  {Bh}  are  uniformly  L2  norm 
equivalent  if  there  is  a  constant  c  >  0  independent  of  h  such  that 

iAhBh\^c    ,     \BhAh\^c    .  (1.15) 

We  say  that  A^^  and  Bh^  are  uniformly  L2  norm  equivalent  if  there  is  a  constant  c  >  0 
independent  of  h  such  that 

\Bh^AH\L,^c   ,     [Ah^BhV^^c    .  (1.16) 

One  sees  immediately  that  (1.15)  is  equivalent  to  (1.11a)  and  (1.16)  is  equivalent  to  (1.12a). 

The  concepts  of  spectral  equivalence  and  L2  norm  equivalence  are  different  even  for 
positive  symmetric  operators.  For  compact  operators,  norm  equivalence  implies  spectral 
equivalence.  The  converse  is  not  true  [FMP].  In  Section  4,  Example  6  we  give  a  simple 
example  of  two  self-adjoint  positive  definite  ordinary  differential  operators  M  and  N  which 
are  spectrally  equivalent  while  M~^  and  yv~^  are  not  L2  norm  equivalent  and  M  and  N  are 
not  Li  norm  equivalent.  Similarly,  if  we  let  {Mh}  and  {A^^}  denote  the  usual  finite-difference 
approximations  to  M  and  A^  respectively,  {Mh}  and  {A^^}  are  uniformly  spectrally  equivalent 
while  neither  {Mh^}  ,  {l^hH  nor  {Mh}  {Nh}  are  uniformly  L2  norm  equivalent. 
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The  Hi  norm  equivalence  A  ^  and  B  ^  is  closely  related  to  spectral  equivalence  of  the 
forms  <Au,v>L  and  <Bu,v>l  when  A  and  B  are  positive  definite  and  self-adjoint  in 
Liiil).  For  general  operators  we  say  that  A"^  and  B~^  are  Hi  norm  equivalent  on  12(0)  if 
there  exists  a  dense  set  D  C  Lii^)  and  constants  o  <  ci  ^  C2  such  that  for  every 
f  ^  OJ  ^  D  we  have 

Although  definition  (1.17)  will  suffice  for  our  purposes,  we  choose  to  reformulate  it  as  a 
stronger  statement  about  the  weak  forms  of  the  operators  A  and  B .  In  Section  2  Aj^  and  S„, 
are  defined  on  HiiClJoiA))  Q  HiiO.)  and  Hi(il,ro{B))  C  //(H)  respectively.  We  say  that 
Aw^  is  //i  norm  equivalent  to  B^^  if  there  exists  Di  C  W^iCn.FoCS)),  Di  dense  in 
HiinjoiB)),  D2  C  //iCnjoCA)),  D2  dense  in  ^iCnjoCA)),  and  constants  ci,  C2>  0  such 
that  for  every  v#0,v€Diwe  have 

Uv^Xvk  :s  cilvk  (1.18a) 

and  for  every  u  #  0,  m  €  D2  we  have 

l5v7^A^«k<  C2lwk   •  (1.18b) 

As  stated  earlier,  bounds  of  the  form  (1.7)  are  not  necessary  for  the  Hi  results.  Rather,  we 
restrict  ourselves  to  operators  for  which  a  special  relationship  exists  between  ao  and  ai  in 
(1.1c).   This  will  be  explained  in  Section  2. 

In  Section  2  we  recall  some  basic  facts  about  elliptic  operators  and  their  adjoints.  We 
also  discuss  the  "weak"  forms  of  elliptic  operators.  In  Section  3  we  discuss  norm 
equivalence  of  elliptic  operators.   The  basic  results  are 

Theorem  3.1.  Let  A  and  B  be  two  invertible  uniformly  elliptic  operators  of  the  type 
described  by  (1.1).  Assume  the  domain  Xl  is  smooth  and  all  invertible  uniformly  elliptic 
operators  of  the  form  (1.1)  whose  coefficients  are  smooth  and  whose  boundary  conditions 
are  based  on  {ro(A),ri(A)}  or  {ro(S).ri(5)}  yield  estimates  of  the  form  (1.7),  (1.9).  (For 
example,  assume  (i)  either  ToiA)  =  0  or  TiiA)  =  0,  and  (ii)  either  ro(5)  =  0  or  TiiB)  = 
0).   Then, 

I).     A"^  and  B~^  are  Z.2  norm  equivalent  on  L2(fi)  if  and  only  if  D(A*)  =  D(B*), 

II).    A  and  B  are  L2  norm  equivalent  on  LziH)  if  and  only  if  D(A)  =  D(B).    ■ 

Theorem  3.2.  Let  A  and  B  be  two  invertible  uniformly  elliptic  operators  of  Class  A^.  Then, 
A^^  and  fi;;^  are  Hi  norm  equivalent  on  LiiO,)  if  and  only  if  ro(A)  =  ro(5);  that  is,  there  is 
a  constant  K^iA-.B)  such  that 

Ujifl,vvk  ^  K2(A:B)  Ivk 
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for  every  v  €  Hii^XoiB))  and 

for  every  v  €  //iCn.FoCA))  if  and  only  if  TqCS)  =  Tq{A).    ■ 

In  Section  4  we  discuss  some  examples  which  illustrate  the  results  of  Section  3  and  pro- 
vide insight  into  the  results  of  Section  5.  In  Section  5  we  discuss  the  uniform  Li  norm 
equivalence  of  {A^"^}  and  {Bh'^},  the  uniform  L2  norm  equivalence  of  {A;,}  and  [Bh)  and  the 
uniform  Hi  norm  equivalence  of  {Ah'^]  and  [Bh'^}.  In  all  cases,  such  uniform  norm 
equivalence  of  the  discrete  operators  can  occur  only  if  the  corresponding  norm  equivalence 
holds  for  the  operators  A  and  S  or  A"^  and  5"^  (see  [FMP]).  Moreover,  under  reasonable 
assumptions  (including  optimal  order  convergence)  the  L2  norm  equivalence  of  A~^  and  B~^ 
implies  the  uniform  L2  norm  equivalence  of  {Ah^}  and  {5a" ^}  and  the  Li  norm  equivalence  of 
A  and  B  implies  the  uniform  L2  norm  equivalence  of  [Ah)  and  {Bh}.  The  first  of  these 
results  was  stated  in  [BP].  Finally,  for  a  particular  class  of  finite-element  methods,  the  Hi 
norm  equivalence  of  A"^  and  B~^  implies  the  uniform  H\  norm  equivalence  of  [Ah^  and 
[Bh'}. 

We  close  this  introduction  with  the  observation  that  one  would  usually  combine  precon- 
ditioning with  a  conjugate  gradient  type  iteration.  Bounds  on  Ci^{Bh^Ah),  CiJ^AhBh^) 
merely  yield  upper  bounds  for  the  rate  of  convergence  of  conjugate  gradient  methods.  A 
thorough  analysis  of  CG  methods  requires  more  information  on  the  distribution  of  either  the 
eigenvalues  or  the  singular  values  of  S/T^A/,  depending  on  the  implementation.  We  discuss 
these  questions,  together  with  some  detailed  computational  experiments  in  another  paper. 

2.   Background 

The  discussion  in  this  section  recalls  some  basic  facts  about  elliptic  differential  opera- 
tors. This  discussion  clarifies  the  facts  concerning  "natural  boundary  conditions,"  the  role 
of  tangential  derivatives  and  the  relationship  to  the  boundary  conditions  of  A*.  In  addition 
we  discuss  some  of  the  basic  facts  concerning  A„,,  the  "weak"  or  Hi  form  of  the  operator  A. 
This  discussion  depends  on  the  "Green's  identities"  for  the  differential  operator  A. 

With  the  operator  A  described  by  (1.1a)  and  the  boundary  conditions  (1.1b),  (1.1c)  we 
also  have  the  operator  A  given  by  (1.1a)  without  any  boundary  conditions. 

Let  injc,ny)  =  n{x,y)  be  the  unit  outward  normal  to  dil  at  a  point  {x,y)  €  30.    Then 

-^  =  Uj,n^  +  UyTiy    ,  (2.1a) 

-^  =  -Uj^riy  +  Uyrtj,   ,  (2.1b) 

The  conormal  derivative  is  given  by 

Ultracomputer  Note  143  Page  7 


"37"  •  ^  t'^ii"'  "*"  ^^"yl^^c  +  [anux  +  a22W>]«>    •  (2.1c) 

This  conormal  derivative  represents  the  "flux"  associated  with  the  operator  A. 

We  assume  that  the  divergence  theorem  holds  in  fl  and  immediately  obtain  the  first 
"Green's  Identity":   let  u  and  v  be  smooth  functions  defined  on  ft.   Then 

iAu,v)  =  a(M,v)  -  /  V  -^  da  (2.2a) 

where 

a(M,v)  =  aiiu,v)  +   ff  v[aiu^  +  azUy  +  <3o"]^^^y  (2.2b) 

and  ai{u,v),  the  bilinear  form  associated  with  the  "leading  part"  of  A  is  given  by 

aL(".v)  =  ^  [anUxVx  +  ai2(«xVy  +  UyV^,)  +  a22"yVy}Jj:Jy    .  (2.2c) 

One  more  application  of  the  divergence  theorem  yields 

(,Au,v)  =  aL(u,v)  -  (w.[(aiv)^  +  ia2v)y])  +  (aou,v) 

-  J    [v  -^^  -  (ain^:  +  a2«>')"v]Jcr    . 

Let  A*  denote  the  (formal)  adjoint  of  A.  That  is,  there  are  no  boundary  conditions  associ- 
ated with  A*  and 

-^'v  =    -  {(aiiV;t)^  +  (ai2V;y);t  +  (flnvjy  +  (a22Vy)>} 

(2.4) 
-  (aiv).t  -  (a2v);y  +  aov    . 

The  conormal  derivative  depends  only  on  the  leading  part.  Thus,  the  conormal  associated 
with  A*  is  the  same  as  the  conormal  derivative  associated  with  A.  One  more  application  of 
the  divergence  theorem  yields  "Green's  Second  Identity." 


(A 


M.v)  -  (w,A'v)  =    f 


bv  du 

U    -T V    -: 

dVA  dVA 


d(T 

(2.5) 


+  Jiairix  +  a2ny)uvdc 


The  operator  A*,  the  L2  adjoint  of  A,  is  obtained  by  using  the  boundary  conditions  (1.1b), 
(1.1c)  and,  if  necessary,  integration  by  parts  along  Fj  (converting  tangential  derivatives  of  u 
into  tangential  derivatives  of  v)  to  determine  appropriate  boundary  conditions  for  A*  so  that 

(Au,v)  =  iu,A*v)    .  (2.6) 

That  is.  A*  is  given  by 

A*v  =  A*v     in  ft  (2.7a) 
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and,  on  dD,,  there  are  prescribed  boundary  conditions  of  the  form 

V  =  0    on  To   .  (2.7b) 

1^  =  a^v  +  ai*  j^     on  Ti    ,  (2.7c) 

where,  ao  and  aj  have  been  chosen  so  that  (2.6)  holds. 

In  order  to  find  ao,aJ  it  is  convenient  to  rewrite  -r—  as 

3«    _  w   du    ,   J,   du 


dv^  dv  dcT     ' 

where 

N  =  a\in}  +  2ai2«;c«y  +  022iy  —  ^(A)    , 

T  =  anin}  -  nj)  +  (^22  "  an)nxny    . 

Observe  that  -j-^  is  always  a  nontangential  derivative  (by  (1.2)  N  ^  X(A)  >  0).    When  the 

leading  part  of  A  is  a  diffusion  operator  with  an  =  aii  and  (2i2  =  0  the  conormal  derivative 
is  a  positive  multiple  of  the  normal  derivative.  Moreover,  unless  one  imposes  artificial  rela- 
tionships between  (n^./iy)  and  the  coefficients  an,  a\2,  022.  this  is  the  only  class  of  operators 
for  which  this  is  so. 

In  many  problems  the  boundary  conditions  associated  with  the  operator  A  are  of  the 
special  form 

-^  =  a(o-)M     on  Ti   .  (2.8) 

Indeed,  the  "natural"  boundary  condition  associated  with  A  are  given  by 

du 


3va 


=  0    on  Ti   .  (2.9) 


When  the  boundary  conditions  for  A  are  of  the  form  (2.8)  we  may  use  (2.5)  to  see  that  the 
boundary  conditions  for  A*  are  given  by 

-^  =  [a  -  {airij,  +  a2«j,)]v:  =  a*v    .  (2.10) 

These  formulae  show  why  our  discussion  must  generally  include  tangential  derivatives  and 
how  the  coefficients  ai,a2  effect  the  boundary  conditions  of  A*. 

In  general,  A  is  self-adjoint  (A  =  A*)  if  and  only  if 

(i)     fli  =  a2  =  0 

(ii)     -^  =  au  on  Fi. 
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As  we  have  written  the  operator  A  in  (1.1a)  one  expects  to  deal  with  a  theory  involving 
second  derivatives,  i.e.,  an  H2  theory.  However,  it  is  well  known  that  there  is  an  H\  theory 
based  on  (2.2)  (cf  [Ci]).  Indeed,  this  theory  is  the  basis  of  much  of  the  finite-element 
method.  While  there  are  some  discussions  of  the  finite-element  method  to  general  boundary 
conditions  (see  [BP])  much  of  that  theory  is  limited  to  the  case  where  the  boundary  condi- 
tions are  given  by  (1.1b)  and  (2.8).  We  will  denote  such  operators  as  operators  of  "Class 
N"  to  indicate  that  the  boundary  conditions  involve  only  the  conormal  derivative.  Our  dis- 
cussion of  the  H\  theory  and  H\  norm  equivalence  will  be  limited  to  this  class  of  problems. 

Such  problems  can  be  reformulated  as  follows.    Let  Fq  =  Tq(,A)  and  let 

//i(n,ro):  =  {(}>  ^  //i(n):7o4)  =  O     on  To}  (2.11) 

where  70  denotes  the  "trace"  operator  (see  [Gr]).    The  function  which  satisfies  (1.1b),  (2.8) 
and 

Au  =  f    inn  (2.12) 

also  satisfies 

a(M,v)  -      f    CLUvda  =  (/,v)    ,     Vv  €  HiiClJoiA))    .  .^  ,-,. 

Equation  (2.13)  represents  the  "weak"  form  of  the  boundary-value  problem  (2.12),  (1.1b), 
(2.8). 

Let  us  reinterpret  (2.13).  Given  a  boundary  value  problem  involving  an  operator  of 
Class  N  we  consider  the  bilinear  form  d{u,v)  on  HiiD-Xo)  ^  HxinXo)  defined  by 

d(u,v)  =  a(u,v)  -   f  auvda    .  ^2  14) 

For  each  fixed  u  €  Hii^.To)  the  value  d(u,  •)  is  a  bounded  linear  functional  on  Hi(D,,ro)- 
The  weak  form  A^  of  A  is  the  mapping  taking  u  into  d(u,  •).   That  is, 

(A^u)(v)  =  d(u,v)    .  (2.15) 

Thus,  (2.13)  is  the  statement  that  (A^^u)  ^  L2i^)  in  the  sense  that  A^u  may  be  represented 
by /as 

(Aw")(v)  =  (f,v)  (2.16) 

We  have  given  this  careful  and  explicit  review  of  these  ideas  to  emphasize  the  fact  that 
in  this  Hi  formulation  of  boundary  value  problems;  if  A  and  B  are  two  uniformly  elliptic 
operators  of  Class  A^  and  ro(A)  =  ro(5),  then  A„,  and  5„,  have  the  same  domain,  Hi(D,,ro). 

When  working  with  d(u,v)  one  must  be  able  to  bound  boundary  integrals.  The  basic 
estimate  is:  there  is  a  constant  Lq  =  Lo(D,),  depending  only  on  H,  such  that,  for  every 
V  €  Hiin), 
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f  v2ja  <  Lo  [iVvL^  •  IvL,  +  |v|^^]    ,  (2.17) 

(see  [Gr]).  In  fact,  it  is  this  estimate  which  enables  one  to  assert  that  A„,m  is  a  bounded 
hnear  functional  on  Hi(il,To).  That  is,  there  is  a  constant  Li(A)  >  0,  such  that,  for  all 
M,v  €  f/i(n,ro)  we  have 

|fl(M.v)|<Li(A)|M  \h^-  \v\h,   .  (2.18) 

If  diu,v)  is  symmetric,  i.e.,  d(u,v)  =  a(v,M),  (ai  =  a2  =  0)  we  say  Ay^  is  positive 
definite  if  there  is  a  constant  L2  =  LaCA)  >  0  such  that 

L2(A)|v|^^  <  a(v,v)    ,     V  V  €  //i(n.ro)    .  (2.19) 

Observe,  that  if  A^  and  S^,  are  two  such  positive  definite  operators  with 
To(A)  =  ro(5),  then,  not  only  do  A^v  and  S„,  have  the  same  domain  but  they  are  spectrally 
equivalent  because 


3.   The  Differential  Operators 

In  this  section  we  discuss  the  basic  facts  about  L2  norm  equivalence  of  elliptic  differen- 
tial operators  as  well  as  some  results  on  the  Hi  norm  equivalence  of  A^^  and  B^^.  The 
basic  results  are  Theorem  3.1  and  Theorem  3.2. 

Lemma  3.1.  Let  A  and  B  be  two  invertible  uniformly  elliptic  operators  of  the  form 
(l.la,b,c)  such  that  (1.7)  and  (1.9)  hold.  Let  KM)  Ki(B),  KjiA),  KjiB)  be  the  constants 
associated  with  the  bounds  (1.7),  (L9).   Assume 

D{A)  =  D{B)    .  (3.1) 

Then,  AB~^  and  BA~^  are  defined  on  all  of  L2(ft)  and 

\AB-\^  K2{A)Kx{B)    ,  (3.2a) 

\BA-\^  Ki{A)K2m    .  (3.2b) 

Thus,  A  and  B  are  L2  norm  equivalent. 

Proof.   The  proof  is  immediate.  ■ 

Lemma   3.2.     Let  A    and   B    be   two   invertible   uniformly   elliptic  operators   of  the   form 

(Lla,b,c)  such  that  (1.7)  and  (1.9)  hold  for  A*  and  B' .    Let  Ki(A'),  K2(A*).  Ki{B*),  K2(B') 

be  the  constants  associated  with  the  bounds  (1.7),  (1.9).   Assume 

DiA')  =  D(B*)  (3.3) 
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Let 

D:  =  D(A)  n  D(B)    .  (3.4) 

Then,  D  is  dense  in  Lzi^)  and  the  estimates  (1.13)  hold  in  the  form 

]B-W\l,^  Ki(B')K2iA')\f\L^   ,    fiD{A)    ,  (3.5a) 

[a-W\l,^  Ki(A')K2(B*)\f\L,   ,    f^DiB)    .  (3.5b) 

Hence,  A~^  and  S~^  are  L2  norm  equivalent. 
Proof.   The  set  D  is  dense  because 

c^{n)  CD. 

Consider 


Hence,  (3.5a)  follows  from  Lemma  3.1.   A  similar  argument  yields  (3.5b).   ■ 

We  now  turn  to  the  construction  of  special  functions  w.^,  (j),,  which  enable  us  to  prove 
the  necessity  of  (3.1)  and  (3.3).  These  functions  will  be  bounded  in  L2  but  become  very 
large  on  a  small  segment  of  the  boundary  where  the  boundary  conditions  of  A  and  B 
disagree.  For  the  next  two  lemmas  we  assume  that  ft  has  a  particular  shape  (see  Fig.  1). 
Condition  ft.  ft  is  contained  in  the  right-half-plane  x  >  0  and  is  convex.  Moreover,  there  is 
an  interval 

T2:  =  {(0,y):[y|  <  8}C  aft 

and  the  entire  rectangle 

Ri-  =  {(x,yy.    (yl  ^  y    .     0  <  ;c  <  1} 

is  contained  entirely  in  ft  with 

dRi  n  aft  c  r2  .  ■ 
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Let  viy)  be  a  fixed  function  which  satisfies 

V  €  C^C-oo^)  (3.6a) 

0  <  v(y)  <  1  (3.6b) 

v(y)  =  1    .     \y\^J  (3.6c) 

v(y)  =  0   ,    lyl>  Y  .  (3.6d) 

Lemma  3.3.   Let  S  be  an  invertible  uniformly  elliptic  operator  given  by 

Bu   =     -   {(bnUx)x   +   ibi2Uy)x   +   (fel2«:.)y   +   {b22Uy)y\ 

(3.7a) 

+     &!«;(     +     b2Uy     +     fco" 

M  =  0    on  ro(fi)  (3.7b) 

-f^  =  Po"  +  3i  -f^    on  ri(5)    .  (3.7c) 

Suppose  n  satisfies  Condition  fl  and 

FaCriCS)   . 
Then,  for  every  ti,  0  <  nri  ^  1,  there  are  two  functions  w,,,  cf),,  €  C^CO)  which  satisfy 

support  Uy^C  Ry^    ,    support  ^y^C  Ry^    ,  (3.8a) 

where  ^^i  is  the  rectangle 

R^.  =  {{x,y):   tyl^y   .    0  ^  a:  <  ti}   .  (3.8b) 

The  function  u,,  satisfies 

ur,  €  D(5)    ,  (3.9a) 

M^(0,y)  =  ^l-^'(y)    ;     (0,y)  €  ail   ,  (3.9b) 

I"tiIlj^c    ,     |«^U,  <cTi-l    ,     lu^l//^  <  cTi-2    ,  (3.9c) 

for  some  constant  c  >  0,  independent  of  t|.   The  function  cj),,  satisfies 

<})^  =  0     on  an    ,  (3.10a) 

■^  {Q,y)  =  Ti-3'2  ^(y)  .    (o,>-)  ^  an  ,  (3.iob) 
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l4>r,lL,  =S  C,  |4>^I//,  <  CT]-!     ,      Ic})^!//^  <  CTl-2     ,  (3.10c) 

for  come  constant  c  independent  of  t].    Proof.   Define 


8{x:r\)  = 
Observe  that  g(,x:-f\)  i  C\0,'^)  and 


(l_^)3(i  +  ^)    ,   O^x^^    .  ^3^^^ 


0     ,  Tl    <   X 


g(0:Ti)  =  1    .     -^^(0:^1)  =  0    .  (3.12) 

Set 

w„U,y)  =  ^-''{viy)  -  x[Po(y)v(y)  +  ^i{y)v'{y)]}g{x:y\)    ,  (3.13a) 

^,^iix,y)  =  -r\-^'^v{y)xg{x:-x\)    .  (3.13b) 

The  desired  estimates  all  follow  from  elementary  computations.    ■ 

Lemma  3.4.    Let  fi  and  5  be  as  in  Lemma  3.3.    Let  A  be  an  invertible  uniformly  elliptic 

operator  which  differs  from  B  only  in  the  boundary  conditions.   That  is, 

A  =  B    ,    D{A)  i=  D{B)    .  (3.14a) 

Specifically,  let  A  satisfy  the  boundary  conditions 

M  =  0    ,     on  Tq{A)    ,  (3.14b) 

1^  =  aoM  +  ai  -f^    ,     on  Y^{A)    .  (3.14c) 

We  consider  several  cases. 
Case  1.    T2  C  ri(A),but 

lao  -  Pol  +  Ui-Pil  >  0    on  r2   .  (3.15) 

Let  Mt,  be  the  function  given  by  (3.13a).   Then,  there  is  a  constant  ci  >  0  such  that 

\A-^Bu^\i^^  ciy\-%iy^\i^   .  (3.16) 

Case  2.    Tj  C  Yq{A). 

Then  (3.16)  holds.    Moreover,  let  (j).^  be  the  function  given  by  (3.13b).   Then 

<t)„  ^  D{A) 
and  there  is  a  constant  C2  >  0  such  that 

]B-^A<^r,\L,^C2-^-^'^\^X   ■  (3.17) 


Ultracomputer  Note  143  Page  14 


Proof.   Let 

A-i5u^  =  Af^    ,  (3.18a) 

w^  =  u^  -  v|j.p    .  (3.18b) 

Then 

and  w.p  satisfies 

Aw^  =  5wt,  =  0    .  (3.19a) 

and  since  i}*^,  =  (w.^  —  m,,)  ^  D(A),  the  boundary  conditions  of  A  applied  to  w.^  is  equal  to 
the  boundary  conditions  of  A  applied  to  u-^. 

This  yields 

Tv^  =  0    on  To(A)    ,  (3.19b) 

^  =  aow^  +  "1  -^     on  ri(A)/r2   .  (3.19c) 

Case  1 .   On  r2  we  have 

(Wt,)x  +  o.QW^  +  ai(w.p)y  =  (m.^);c  +  aoMi,  +  ai(M.p)y  =  •ti"^[.s(}')]  (3.20a) 

where 

s{y)  =  (ao  -  Po)v(y)  +  (ai  -  3i)v'(y)    .  (3.20b) 

Observe  that  v(y)  can  be  chosen  so  that  s{y)  is  not  identically  zero. 
Case  2.   On  r2  we  have 

Wt,  =  Mt,  =  Ti"^'    .  (3.21) 

In  both  cases  the  fact  that  A  is  invertible  means  that  w.^  is  uniquely  determined  by  (3.19), 
(3.20),  and  (3.21).    However,  a  simple  computation  shows  that  ti'^Hvi  is  a  solution.    Hence 

and  the  estimates  (3.16)  follows  from  (3.9c)  and  the  triangle  inequality. 

It  is  clear  from  (3.10a)  that  in  Case  2,  4)^,  ^  D(A).    In  order  to  prove  (3.17),  let 

fl-U(})^  =  vj;^    .  (3.22a) 

w,!  =  c|)^  -  »|/^    .  (3.22b) 

Then 

Aw^  =  Bw^  =  0    .  (3.23) 
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Now,  w^  satisfies  the  boundary  conditions  of  B  on  diyr2,  and  on  r2  we  have 

iw^)x  +  PoH-^  +  Pi(w^)^  =  i<^^),  +  po4>^  +  pi(4>^),  =  Ti-3/2^(y)    .  (3.24) 

Once  more,  the  invertibility  of  B  implies  w^  is  the  unique  solution  of  (3.23)  and  (3.24). 
Again 


w. 


=    T\     ^^^1 


and  we  obtain  (3.17).    ■ 

Theorem  3.1.  Let  A  and  B  be  two  invertible  uniformly  elliptic  operators  of  the  type 
described  by  (1.1).  Assume  the  domain  fl  is  smooth  and  all  invertible  uniformly  elliptic 
operators  of  the  form  (1.1)  whose  coefficients  are  smooth  and  whose  boundary  conditions 
are  based  on  these  partitions  yield  estimates  of  the  form  (1.7),  (1.9).  (For  example,  assume 
(i)  either  ro(A)  =    0  or  ri(A)  =    0,  and  (ii)  either  ro(5)  =    0  or  Tx{B)  =    0).   Then 

I)  A"^  and  B~^  are  L2  norm  equivalent  if  and  only  if  D{A*)  =  D(B'), 

II)  A  and  B  are  L2  norm  equivalent  if  and  only  if  D{A)  =  D{B). 

Proof.  Since  the  "if"  parts  of  these  statements  are  proven  in  Lemma  3.2  and  Lemma  3.1, 
respectively,  we  need  only  prove  the  "only  if"  part.    We  start  with  statement  (I). 

Because  L2  norm  equivalence  is  transitive  we  use  Lemma  3.2  to  conclude  that  we  need 
only  consider  the  case  where 

A  =  B    ,  (3.25a) 

and 

Df^A*)  #  D(5*)    .  (3.25b) 

However,  (3.25a)  and  (3.25b)  imply  that 

D{A)  #  D(5)    .  (3.26) 

Let  Fs  C  dfl  be  a  portion  or  the  boundary  on  which  the  boundary  conditions  are  dif- 
ferent. By  a  smooth  mapping  we  may  map  ft  onto  a  domain  ft'  which  satisfies  Condition  ft 
and  r2  C  Fs,  the  image  of  F3.  Using  Lemma  3.4  we  see  that  A'^B  and  S"U  are 
unbounded  on  their  domains.  Hence  they  cannot  be  bounded  on  any  dense  set.  Thus,  we 
have  proven  statement  (I). 

Consider  statement  (II).  Once  more  it  suffices  to  consider  the  case  where  (3.25a)  and 
(3.25b)  hold.  The  proof  proceeds  by  contraposition.  Assume  there  is  a  dense  set  D  for 
which  (1.13)  holds.    For  every/  €  D  we  have 

Ufi-i/L,  =  sup    l-^^^^l^    ,     V  €  L2,v  ^  o|   .  (3.27a) 

Since  D{A*)  is  dense  in  L2(^)  we  may  restrict  the  function  v  of  (3.28a)  to  be  in  D(A*);  that 
is, 
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US-I/L,  =  sup    \^^^;\f^    .     v-  €  DiA*),v  #  ol   .  (3.27b) 

Let  0  <  €  <  1  be  given.    Since  statement  (I)  holds  and  (3.25b)  holds,  there  is  a  v^  €  D(A*) 
such  that 

\v^L,  =  1    .  (3.28a) 

|(S*)-U*vJ^,  >  I    .  (3.28b) 

Since  D  is  dense,  there  is  an  /^  €  D  such  that 

|(S')-U'v,  -  /4  ^  62  |(S*)-U'vJ^^   .  (3.29) 

Using  (3.28b)  and  (3.29)  we  have 

Ufi-y.L,  2=  I(fi*)-U*veb  [1  -  «^^    •  (3.31) 

Using  the  triangle  inequality  we  have 

KI<(1  +  e2)  |(5*)-iA*vJ   . 

Hence,  (3.31)  yields 

Thus,  (1.13)  cannot  hold  and  the  Theorem  is  proven.   ■ 
Remark.     Observe   that  in   the   proof  of  (I)   we  exhibit  a   sequence  on   which   (S"U)    is 
unbounded  while  in  the  proof  of  (II)  we  prove  that 

F:  =  {f^  Lr.B-^f  ^  DiA)} 

is  QQl  dense.    Both  AB~^  and  BA~^  are  bounded  on  F.    To  see  this  we  merely  observe  that 
inequalities  (1.7)  and  (1.9)  yield 

l4fi-yit,  s  K2(A)Ki(B)\f\L^   ,    f  ^F    . 

Thus,  AB~^  is  bounded  wherever  it  is  defined,  but  it  is  not  defined  on  a  dense  set  unless 
D(A)  =  D{B).    ■ 

We  now  turn  to  the  H\  norm  equivalence  of  A^^  and  B^^  where  both  A  and  B  are  of 
Class  A^.    We  begin  with  a  basic  lemma  on  the  H\  closure  of  D{A).    We  assume  A^u  =  f 
possesses     a     unique     solution     u     for     all     /  €  LiiO).      In     this     case     we     define 
D{Ay  =  {u  =  A-^f;f  €  12(0.)}.   It  may  happen  that  D(A) 
C  Hii^)-    However  we  always  have  (i)  D(A)  C  HiiCLXoiA)),  and  (ii)  there  is  a  8  >  0  and 
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D(A)  C//i  +  8(fl). 

Lemma  3.5.   Let  V  denote  the  Hi  closure  of  D(,A). 

Then 

V  =  Hi(il,To(A))    .  (3.33) 
Proof.   Since  D(A)  C  Hii^JoiA))  (and  since  f/i(fl,ro(A))  is  closed  in  Hi)  we  have 

V  C  HiinjoiA))    . 

Suppose  that  there  is  a  (J)  ^  HiinXoi^))  and  ^  ^  V. 

Consider  the  invertible  uniformly  elliptic  operator  B  of  Class  A^  given  by 

Bu  =    -  [(aii«Jj  +  iai2Uy)x  +  (auujy  +  (a22«y)>}  +  Ku  (3.34a) 

u  =  0  on     ro(A)    ,  (3.34b) 

du 


dVg 


=  au   on  TiiA)    ,  (3.34c) 


where   a   is   the   coefficient   in   the   definition   of  the   boundary   conditions   for  A.     Since 

du  du  , 

"3 —  =  -T —  we  have 

D(A)  =  D{B)    .  (3.35) 

The  constant  K  is  chosen  positive  and  large  so  as  to  ensure  that  B  is  invertible  and  that  the 
bilinear  form 

b{u,v)  =  J  J  (aiiu^v^  +  auiUjcVy  +  UyV^)  +  a22UyVy)dxdy 

.  (3.36) 

+  K{u,v)  —     J     auvdd 
r,(A) 

is  an  //i(ft,ro(A))  inner-product.    That  is,  there  are  constants  0  <  ci  s  C2  such  that 

ci  \u\h^  ^  Hu,u)'^  ^  C2\u\h^   ,     V  w  €  HiinjoiA))    . 

The  estimate  (2.17)  can  be  used  to  show  that  such  constants  exist  for  all  K  sufficiently  large. 
Thus,  b  induces  a  norm  that  is  equivalent  to  the  Hi  norm.  We  may  consider  V  to  be  the  b 
closure  of  D(A). 

Since  4>  ?  V  there  is  a  vj;  €  HiiClXoiA))  which  is  "orthogonal"  to  V  in  the  bi  ■  ,  ■  ) 
inner  product.   That  is,  for  all  m  €  D(A) 

b{u,^)  =  0    . 

However,  for  each  /  €  L2(^)  we  may  solve 

Bu  =  f   .  (3.37) 
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Clearly,  the  solution  u  ^  D{B)  =  D{A).   But,  the  discussion  of  Section  2  shows  that 

0  =  b(uM  =  (f,^\,)  (3.38) 

Hence,  ^  €  //i(fi,ro(A))  is  L2  orthogonal  to  Lii^),  which  implies 

4/  =  0 

and  the  lemma  is  proven.   ■ 

The  main  result  which  follows,  namely.  Theorem  3.2,  can  now  be  proven  by  a  rather 
general  argument.  However,  we  proceed  with  the  following  more  detailed  discussion  as  it 
immediately  generalizes  to  the  discrete  case. 

Lemma  3.6.    Let  A  and  B  be  two  invertible  uniformly  elliptic  operators  of  Class  A^  which 
satisfy 

Au  =  Bu  +  CoU  (3.39a) 

where  Co  t^  0  is  a  constant,  and 

D{A)  =  DiB)    .  (3.39b) 

Let  u,v  €  DiA)  =  D(S)  satisfy 

Ah.m  =  B^v  (3.40) 

Then,  there  is  a  constant  K^CA-.B)  such  that 

l"k  -  ■^3(A:5)  \v\h^   ,  (3.41a) 

\v\h^^  Ki{A:B)  \u\h^   .  (3.41b) 

Proof.   Let 

w  =  u  —  V 
then 

Aw  =  Au  —  [Bv  4-  C(,v]  =    —  CoV    . 
That  is 

Aw  =    —  CqV    . 

Since  v  ^  D{A),v  €  L2(fi)  and  the  invertibility  of  A  on  Li^D.)  implies  that  A„,  is  invertible  on 
Lii^)  and 

Vv\h,  ^  Ki[:o\  \v\l^  <  Ki\co\  Iv|^^ 

for  some  appropriate  constant  Ki  >  0.    Similarly 

Bw  =  [Au  —  cqu]  —  Bv  =    —  cqu 

and  for  some  K2>  0 
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The  estimates  (3.41a),  (3.41b)  now  follow  from  the  triangle  inequality.     ■ 
Corollary.   Let  A„,  and  B^  be  as  in  Lemma  3.6.    If  u,v  ^  HiiilXoi^))  satisfy 

Then  (3.41a)  and  (3.41b)  hold. 

Proof.  From  Lemma  3.5  we  know  that  D(A)  =  D(B)  is  dense  in  Hi  (n,ro(A)).  Thus,  there 
exists  a  sequence  u„  €  D{A)  such  that  \u„  —  u\h  ->  0.  Since  B  is  invertible,  there  exists 
v„  €  D(S)  such  that 

A^U„   =  AUr,   =   Bv„   =   SvvVfl 

for  every  n.  Now,  Lemma  3.6  implies  that  |v„  —  v|//^  -  0  for  some  v  €  HiiflXoiA)).  The 
boundedness  and  invertibility  of  A^  and  B^  yield  v  =  v.  Applying  (3.41a,b)  to  «„,  v„  and 
taking  limits  yields  the  result.    ■ 

Theorem  3.2.  Let  A  and  5  be  two  invertible  uniformly  elliptic  operators  of  Class  N .  Then, 
AZ^  and  B^^  are  Hi  norm  equivalent  on  LiiO,)  if  and  only  if  VoiA)  =  FoiB);  that  is,  there  is 
a  constant  K:iiA:B)  such  that 

U>:XvI//^  <  /i:3(A:5)  IvW^  (3.42a) 

for  every  v  €  HiiajoiB))  and 

\B-^A^u\h^^  K3iA:B)\u\H^  (3.42b) 

for  every  v  €  HiinJoCA))  if  and  only  if  ro(5)  =  ro(A). 

Proof.    We  will  prove  (3.42a).    Since  norm  equivalence  is  transitive  we  may  use  Lemma  3.6 

and  assume  that  the  coefficient  oq  of  the  operator  A  satisfies 

aQ>:K  (3.43) 

for  any  suitably  chosen  K  >0  .   We  shall  choose  K  later. 

Assume  A^u  =  B„v  and  To(A)  =  ro(5).    Using  the  representations  (2.14),  (2.15)  we 
have  for  every  <\>  d  Hi  (ClXoiA)) 

a(w,4>)  -    r  0LU<\)d(T  =  b(v,4>)  -  J  ^v^dd    .  ^2  44) 
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Setting  4)  =  u  and  using  the  form  of  A  and  B  this  equation  yields  an  estimate 


+  cil 


+  C2 


r  u^da-  +  J  v^dc 


(3.45) 


where  the  constant  ci,  depends  on  aj,  a2,  i^i,  bz,  bo  but  not  on  ao  while  the  constant  C2 
depends  on  a  and  3- 

Using  the  inequality 


ab^^a^+^b^ 


(3.46) 


several  times  in  (3.45)  we  obtain  the  estimate 

I  X(A)|Vu|^.  +  I  KlA^^^  <  c  [iVvl^.  +  Ivl^;] 


(3.47) 


+  C2 


r  «2jct  +  r  v2j(T 


where  c  is  a  large  constant.    Using  the  estimate  (2.17)  we  may  omit  the  two  boundary 
integrals  and  use  larger  values  of  c  and  ci.   That  is,  we  now  have 


\  X(A)  |VuI^5  ^  T  ^'"'i:  -  ^  ''''^1  +  ^1  '"'^2  ■  '^"'^2 


(3.48) 


where  c  depends  on  ai,  a2.  ^i.  ^2.  ^O  and  Lq  =  LQ{il),  but  c  does  not  depend  on  K.    We 
now  choose 


K  = 


4c? 


1  AMI 

X(A)     '     ^        2      ci 


(3.49) 


and  apply  (3.46)  once  more.   This  yields 

and  (3.42a)  follows  with  K^^A-.B)  =  c7min  {j-  \(A),  -^). 

Now  suppose  ro(A)  n  Yi{B)  #  0.  By  transitivity  we  may  assume  A  =  B.  Let  us  map 
n  onto  a  domain  with  Condition  fi  so  that  r2  C  ro(A)  fl  ro(S).  Consider  the  function 
<?T)  =  Tl<l>-ri  where  (j)^,  is  defined  in  (3.13b).    By  (3.10c)  we  have  \q-^H^  —  c. 
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Following  the  development  in  the  proof  of  Lemma  3.4  yields 
for  some  c  >  0.   This  completes  the  proof.    ■ 


(3.50) 


4.   Examples 

In  this  section  we  look  at  simple  examples  which  illustrate  the  results  of  Section  3  and 
give  insight  into  the  results  of  Section  5.  The  operators  here  are  all  ordinary  differential 
operators  defined  on  LilO,!].  In  each  case  we  carry  out  explicit  calculations,  based  on  the 
appropriate  Green's  function,  to  exhibit  the  relevant  results  of  Section  3. 

Let 


Lu  =    -  u"    , 
Mu  =    -  u"    , 
Nu  =    -  u"    , 
Ru  =  -u"  +  au'    , 
where  a  ^  0.    Notice  that 

R*u  =    -  u"  -  au' 


u(0)  =  «(1)  =  0    , 
u(0)  =  u'il)  =  0    , 
u(0)  =  au(l)  +  «'(1)  =  0 
u(0)  =  u'(l)  =  0    , 

u(0)  =  aud)  +  u'd)  =  0 


(4.1a) 


(4.1e) 


We  will  examine  the  equivalence  of  these  operators  in  the  L2  norm  and  the  equivalence  of 
their  inverses  in  both  the  L2  and  Hi  norms.   The  Green's  function  for  A^  is  given  by 


Gn(x,s)  = 


s(l  - 
xil  - 


1  +  a 
a 


■X) 


1  +  a 


s) 


Q  ■<  s  <  X 
X  <  s  <  \    . 


(4.2) 


The  Green's  functions  for  L  and  M  are  found  by  setting  a  =  <»  and  a  =  0  in  (4.2)  respec- 
tively.  The  Green's  function  for  R  is  given  by 


Gr(x,s)  = 


1 


—  (1  -  e-' 

a  ^ 

(e"  -  l)e- 


1 


0  <  5  <  X 
a:  <  5  <  1 


(4.3) 


The  Green's  function  for  R*  is  found  by  reversing  the  roles  of  ;c  and  5  in  (4.3). 

Example  1).    In  this  example  we  consider  M  and  R.    Since  D{M)  =  D(R),  Theorem  3.1 

implies  that  M  and  R  are  Lj  norm  equivalent.    Since  D{M*)  #  D(R*),  A/"^  is  not  L2  norm 

equivalent  to   /?"^     However,   ToiM)  =  ToiR)    implies   that  A/"^  and   R~^   are   Hi  norm 

equivalent. 

Example  I. a).   L2  norm  equivalence  of  M  and  R.   We  have  D(A/)  =  D(R). 
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A  straightforward  calculation  yields 

1  1 

MR-'^f  =  M  J  Gii(x,s)f(s)ds  =  f  -  ae'"  J  e-'"f(^s)ds 


and 


RM- 


■1/  =  /?!  GM{x,s)f{s)ds  =f  +  a  Sfis)ds 


Since  D(M)  =  D(R),  (4.4a,b)  are  defined  for  all/  ^  L2[0,l].    Write 

MR-^f^f-  ^Ux,s)fis)ds    . 


(4.4a) 


(4.4b) 


(4.5a) 


where 


Then,  we  see 


^l(;c 


_    (       0,  0<5:S;c 


1   1 

r  r  k}{x,s)dxds 


-Vf . 


This  yields 


A  similar  calculation  yields 


lA//?-i/L  ^ 


\RM-^f\L 


'l  +  V  f 


IfL, 


1  + 


vT 


Ifl. 


(4.5b) 


(4.6) 


(4.7) 


(4.8) 


Thus,  as  predicted  by  Theorem  3.1,  we  have  shown  that  M  and  R  are  Lj  norm  equivalent. 

Notice  the  norms  grow  with  the  parameter  a. 

Example  1  .b)  Li  norm  equivalence  of  M~^  and  R~^.   Now  consider  M~^  and  R~^.   We  have 

1 
R-'^Mu  =  jGR(x,s)(-u"is))ds 

J  (4.9a) 

=  u  -  u(0)  +  -J  e-''(e'"  -  l)iau(l)  +  w'(l))  -   f  k2ix,s)uis)ds    , 


where 


k2{x,s)  = 


-ae~°' 


0  ^  s  ^  X 
a(e'"  -  De-"'    ,   x  <  .y  <  1    .' 


(4.9b) 
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For  u  €  D(M)  this  becomes 

1 
R~'^Mu  =  u  +  ae-'ie'"  -  1)m(1)  -   fk2{x,s)u(s)ds    .  (4.10) 

Likewise,  we  have 

M-'^Ru  =  u  -  u(0)  +  x(au(l)  -  m'(1))  -  aju(s)ds    .  (4.11) 

For  u  €  D(R)  this  becomes 

M~^Ru  =  u  +  axM(l)  -  afu(s)ds    .  (4.12) 

Both  (4.10)  and  (4.12)  contain  terms  involving  m(1).    This  quantity  is  not  bounded  on 
D(M)  =  DiR). 

To  see  this  let 

p(z)  =  (4z3  -  3z^)  (4.13) 

The  polynomial  p(z)  has  the  properties 

p(0)  =p'(0)  =  p"(0)  =  0 

p(l)  =  1    ,    p'(l)  =  0 

^max^  \p{z)]  =  1 

If  we  let 


0   ,  0  <  a:  <  1  -  Ti    , 

^-.^  (X^lJ_±JQ.)    ^     i_.^<^<i    ^  (4.14) 


"■nU)  = 

Then,  u^  ^  D(M)  =  D(/?)  and 

K\l,^  1    ,     k(l)l=  y\-'^   . 
Thus,  there  exist  ^i  K2,  tio  >  0  such  that  for  0  <  ti  <  tiq 

\R-'Mu^\L,^Kx'q-'^\u^l^   ,  (4.15a) 

y-^Ru^\L,^K2-n-''\uy,\L,   .  (4.15b) 

Example  l.c).  Hi  norm  equivalence  of  A/"^  and  R~^.  On  the  other  hand,  we  see  that 
ro(A/)  =  FoiR).  Theorem  3.2  predicts  that  \R-^M\h^  and  W~^R\h,  are  both  bounded.  We 
will  use  the  following  lemma. 
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Lemma  4.1.   Let  w  €  //i[0,l].   Then 

|m(1)|<  ^\u\h^   .  (4.16) 


Proof.   We  have 


Thus, 


1 
a(l)  =  u{x)  +  S  u'{x)ds 


|u(l)|<  |«U)|4-   hu'{s)\ds 

^  \ii{x)\  +  lu'li:^   . 

Integrating  both  sides  yields 

1 
I«(l)|<  ]\u{s)\ds  +  \u\ 

Also,  we  note  that  the  integral  operators  in  (4.10)  and  (4.12)  are  bounded  in  the  Hi  norm. 

Taking  the  Hi  norm  of  both  sides  of  (4.10)  and  (4.12)  and  using  Lemma  4.1  establishes  the 

bounds. 

Example  2).    In  this  example  we  consider  A^  and  R.    Since  DiN)  ¥=  DiR),  N  is  not  L2  norm 

equivalent  to  R.   However,  since  D(N*)  =  DiR"),  N~^  is  L2  norm  equivalent  to  R~^.   Also, 

^o(^)  —  ro(^),  so  yv~^  is  Hi  norm  equivalent  to  R~^. 

Example  2. a).   L2  norm  equivalence  ofN  and  R.    Similar  to  (4.4a,b)  we  have 

1 
NR-^f  =  f  -  ae'"  f  e-'"f(s)ds  (4.17a) 

X 

and 

RN-^f  =  f  +  a  Sns)ds  -  j^  J  sf(s)ds    .  (4.17b) 

As  the  bounds  (1.7)  and  (1.9)  imply,  these  operators  are  bounded  whenever  they  apply. 
However,  NR~^  and  R~^N  are  not  defined  on  a  dense  subset  of  L2[0,l].  We  shall  see  in 
Example  5  that  centered  finite  difference  approximations  to  these  operators  yield  discrete 
analogues  to  (4.17a,b)  that  are  not  bounded  as  the  mesh  is  refined. 

Example  2.b).  L2  norm  equivalence  of  N~^  and  R~^.  Since  M  =  N ,  that  is,  M  differs  from  A^ 
only  in  boundary  conditions,  we  have  /?~W  =  R~^M  and  (4.9a,b)  apply.  For  u  €  D{N) 
(4.9a,b)  become 

Ultracomputer  Note  143  Page  25 


1 


R-'^Nu  =  u  -   fk2(x,s)u{s)ds    .  (4.18) 

Similar  to  (4.11)  we  have 

N-'^Ru  =  M  -  (1  -  YT-^x)  u(0)  -    ^  I  ^  xu'il)  -   h2(x,s)u(s)ds        (4.19a) 
where 


ksix.s)  = 
For  u  €  D(R)  this  becomes 


a2 

^  +  ^  2  (4.19b) 

a       r  X  ^  s  <  I    . 


1 


N-'^Ru  =  u  -     ('k3ix,s)uis)ds    .  (4.20) 

Clearly,  R~^N  and  A^~^R  are  bounded  in  Li  norm. 

Example  2.c).  Hi  norm  equivalence  of  N~^  and  R~^.  From  (4.18)  and  (4.20)  we  see  that 
]R~^N\ff^  and  Ia^"^/?!//^  are  bounded.  This  agrees  with  Theorem  3.2  since  ro(A^)  =  FoiR). 
Example  3.  In  this  example  we  consider  L,  M,  and  N .  Notice  that  L  =  M  =  N  and  that 
each  is  a  self-adjoint  positive  definite  operator.  Also  notice  that 
D{L)  #  D(A/)  #  DiN)  i^  D(L)  so  that  no  two  are  L2  norm  equivalent,  nor  are  any  two  of 
their  inverses.  However,  Toil)  ¥"  VoiM)  =  ro(A^).  Thus,  L~^  and  A/"^  are  not  Hi  norm 
equivalent  while  M~^  and  N~^  are  Hi  norm  equivalent.  Finally,  the  bilinear  forms  associ- 
ated with  these  operators,  l(u,v),  m(«,v)  and  niu,v),  are  defined  on 
HiinXod))  i=  //i(n,ro(A/))  =  HiiilJoi^)),  respectively.  Since 

D  =  HiiilXoiL))  n  /fi(fl,ro(M))  is  not  dense  in  either  space,  L  and  M  are  not  spectrally 
equivalent.    On  the  other  hand,  we  will  see  that  m(u,v)  =  n(u,v)  on  //i(n,ro(M)). 
Example  3. a).   Li  norm  and  Hi  norm  equivalence  of  M~^  and  N~^.   Consider 

1 
N-^Mu  =  jGN(x,sX-u"(s))ds 


"  "  "  ^^  "  TTT^^  "^°^  ~  T^T^^''"^^^  ^  "'^^^^ 


(4.21) 


For  M  ^  D{M)  we  have 

/V-^Mu  =  u  -  -j-^A:I^(l)    .  (4.22) 

Recall  that  Gm(x,s)  is  found  by  setting  a  =  0  in  Gn(x,s).   Thus,  A/~W  is  found  by  setting 
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a  =  0  in  (4.21)  to  yield 

M-'^Nu  =  u  -  m(0)  -  xu'{\)    .  (4.23) 

For  u  €  D{N)  we  have  m'(1)  =    —  au{\),  which  yields 

A/-Wm  =  u  +  axu{\)    .  (4.24) 

Both  M~^N  and  N~^M  contain  a  boundary  operator  that  is  unbounded  in  L2  norm  but 
bounded  in  Hi  norm  (see  lemma  4.1).  Thus,  \M~^N\i^  and  |/V~^A/|/,  are  unbounded  while 
]M~^N\h^  and  ]N~^M\h^  are  bounded.  Also  note  that  MN~'^  =  (N~'^M)*  and 
NM-'^  =  (A/-W)*.   Thus,  \MN-\^  and  ]NM-\  are  also  unbounded. 

Example  3.b).  L2  norm  and  Hi  norm  equivalence  of  M~^  and  L~^.  Equation  (4.23)  can  be 
used  to  find  M~^L,  since  N  =  L. 

Here,  however,  u  ^  D{L)  yields 

M-'^Lu  =  u  -  xu'(\)    .  (4.25) 

Similarly,  L~^M  can  be  found  from  (4.21)  by  setting  a  =  °°  to  yield 

L-^Mu  =  u  -  (1  -  x)uiO)  -  xu(l)    . 

For  u  ^  D(M)  this  is 

L-'^Mu  =  u  -  xu(l)    .  (4.26) 

Again,  both  L~^M  and  M~^L  contain  an  unbounded  boundary  term.  Thus, 
\L~^M\l^  =  \ML-\^  and  \M~^L\l^  =  \LM-\^  are  unbounded.  On  the  other  hand.  Lemma  4.1 
shows  that  IL'^mI//^  is  bounded  on  D(M).  However,  W~^L\h^  is  not  bounded  on  D(L).  To 
see  this  let 

q(z)  =  (z^  -  z3)    .  (4.27) 

The  polynomial  q(z)  has  the  properties 

qiO)  =  q'iO)  =  q"{0)  =  0    , 

q(l)  =  0    ,     q'il)  =  1    , 
^max^  \qiz)\  =  ^    . 

imaxy(z)\=}    . 
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If  we  let 


Uj,ix)  = 


0    , 
;c  -  1  +  Ti 

^ 

0  <  X    <    1    -    Tl      , 

1  -    Tl    <  X    <    1      , 


(4-28) 


then  Ml,  ^  D(L)  and 

\uX^f  +  }  .    K(i)l=-n-^ 

Thus,  there  exists  K,t\o  >  0  such  that  for  ti  <  -qo  we  have 


(4.29) 


Example  3.c).    Spectral  equivalence  of  l(u,v),m(u,v)  and  n{u,v).    The  bilinear  forms  associ- 
ated with  the  weak  forms  of  L,  M,  and  A^  are 

1 
l(u,v)  =   f  u'(s)v\s)ds  (4.30a) 


for  u,v  €  Hii[0,l],ro(L))  and 


± 
m(«,v)  =   f  u'(s)v\s)ds 


n(u,v)  = 


f  U'(5)V 


(5)J.s  +  au^d)    ^ 


(4.30b) 


(4.30c) 


for  M.v  €  f/i([0.1],ro(M))  =  //i([0,l].ro(A^)).  Since 

D  =  f/i([0,l],ro(L))  n  f/i([0,l],ro(A/))  is  not  dense  in  either  space,  we  cannot  say  /(m,v)  is 
spectrally  equivalent  to  either  m(u,v)  or  n(u,v).  Example  6.c  will  show  that  the  discrete 
analogues  are  not  uniformly  spectrally  equivalent.    Now  Lemma  4.1  yields 

u(l)2<  2m(u,u)    , 

which  in  turn  yields 


1  <  -^4^  <  1  +  2a 
m{u,u) 


(4.31) 


Thus,  m(u,u)  and  n{u,u)  are  spectrally  equivalent  on  //i([0,l],ro(A/)). 

We  now  turn  to  the  discrete  analogues  of  the  examples  above.    First  we  make  some 
definitions.   Let  n  ^  2  be  a  positive  integer  and  let 

1 


h  = 


n  +  1 


(4.32) 


Define  the  uniform  grid 


Xj  =  jh    ,    ;■  =  0,1,  .  .  .  ,/i  +  1    , 
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on  the  interval  [0,1]  and  consider  the  mesh  functions,  u.  =  {u,}?=o',  that  is,  vectors  of  values 
associated  with  these  points.  We  define  the  discrete  L2  inner  product  between  mesh  func- 
tions 

<U.,y.>h  =  h^UjVj  (4.33a) 

and  discrete  L2  norm 

|uia,A  =  <u.,lL>h    .  (4.33b) 

Observe  that  this  norm  depends  only  upon  the  inner  points  {xj}J=i,  and  differs  from  the  I2 
norm  on  these  points  by  the  factor  h'*.   We  also  define  the  semi-norm 

l«l?,A  =  T  i  (";+l  -  ";)^  (4.34a) 

j  =  0 

and  the  discrete  Hi  norm 

k\lh  =  yi,h  +  kth    .  (4.34b) 

Similar  to  Lemma  4.1  we  have 

Lemma  4.2.    For  any  mesh  function  and  any  0  ^  /  ^  /i  +  1 

Further,  if  uq  =  0,  then 

l"/l  ^  kill,;, 


lulo.A  ^   ^2^  '^'i-'' 


Proof.    The  first  part  is  directly  analogous  to  the  proof  of  Lemma  4.1.    To  see  the  second 
part. 


kP  =  \^  ("y+i  -  w;)P  ^{h^  !)(]-  ^  (Uj+i  -  uj)^) 

j=0  j=0  j=Q 


(4.35) 


^  xi  \ll\lh 

Since  xi  <  \  the  first  bound  is  proven.    Now  using  (4.35)  we  have 

Since  h  =  — -j—r  the  lemma  is  proven.    ■ 

Now  consider  a  centered  finite  difference  approximation  to  the  equation 

Su  =    -  u"  +  au'  =  f   ,     «(0)  =  7;<(1)  +  u'{\)  =  0    .  (4.36) 
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We  have 


Mo  =  0 


ah 


ah 


-  (1  +  ■^)  Uj-i  +  luj  -  (1  -  -^)M;  +  i  =  h^fj    ,  j  =  l,...,n 
-  (1  -  ^)u„  +  (1  +  ^)u„  +  i  =  0    . 


(4.37a) 


(4.37b) 


(4.37c) 


Eliminating  (4.37a)   and  (4.37c)  yields  the  n  x  n   matrix  equation   for  the  interior  mesh 


values,  u.  =  ("1, 


.UnY, 


Shil  =  hH   . 


(4.38) 


where  Sh  is  a  tridiagonal  patrix 
Sh  = 


2.  -  (1  -  ^) 


-  (1  +  ^) 


-  (1  -  ^)  •  •  ■ 

-  (1  +  ■^).  1  +  (3 


(4.39a) 


and 


P  = 


1  - 


(1-^)(1 


ih.^ 


(1  +  ^) 


(4.39b) 


All  of  the  operators  (4.1a-e)  fit  this  pattern.    For  example,  R  has  a  =  a,  7  =  0  which 

yields  p  =  ■^.    On  the  other  hand,  R*  has  a  =    -  a,  y  =  a  which  also  yields  3  =  ■^. 

Thus,  (R*)h  =  (Rh)^  as  one  would  expect.    The  operators  L/,,  Mh  and  A^;,  can  be  similarly 
extracted  from  Si,.   To  find  their  inverses  one  need  only  consider 


.-1  = 


Mh'  = 


1        1         1 

1       2        2 

1       2        3 


1       2 


(4.40) 


and  the  following  Lemma. 

Lemma  4.3.    Suppose  A  =  (B  +  ^mJ)  where  A  and  B  are  both  invertible  n  x  n  matrices. 

Then 


A-i  =  (S-i  -  y^T^ 


where 
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7  =  p/(l  +  P(i^^5-iii)) 
Proo/.   The  proof  is  by  direct  calculation.   ■  Notice  that 

Lh  =  Mh  +  ^e^el  ,       P  =  1    . 


Nh  =  Mh  +  ^e^el   ,       p  = 


ah 


1  + 


ah 


where  ej  =  (0,  .  .  .  ,0,1).    Using  Lemma  4.3,  we  see 

Lh'^  =  A/a"^  -  -Yia:^   ,      y  =  h    , 

ah 


Nf,^  =  Mh  ^  —  yy^^ 


y  = 


1  _i_  oh 

I  +  a  —  -~- 


(4.41a) 
(4.41b) 

(4.42a) 
(4.42b) 


where  yF  =  (1,2,3,  .  .  .  ,n). 
Finally,  we  remark  that 

\Shkh  =  \Sh\i,  ^  \Sh\F  (4.43) 

where  lS;,|;  is  the  Ij  norm  and  IShlp  is  the  Frobenius  norm; 

\A\}=ii  al    . 

With  these  results  we  are  ready  to  examine  the  discrete  analogues  to  the  above  examples. 
Example  4.   This  example  shows  that  Mh  and  Rh  are  uniformly  equivalent  in  the  discrete  L2 
norm  while  Mh^  and  Rh^  are  not,  that  is,  \MhRh\,h  and  I/?a'^a'^Io,a  are  uniformly  bounded 
while  \Mh^Rh\a,h  and  \Rh^MhV),h  are  not.   It  also  shows  that  ]Mh^Rh\i,h  and  \Rh^Mh\i,h  are  uni- 
formly bounded. 

Example  4. a).    Uniform  discrete  Li  norm  equivalence  of  Mh  and  Rh-    Consider  Mh  and  Rh- 
From   (4.36)  and  (4.3^,b)  we  have 


Rh  = 


2    ,  -  (1  -  ^) 
•  •  •    -  (1  +  f") 


-  (1  -  ^)  •  •  • 

-(l  +  ^),(l-f  ^) 


(4.44) 


Multiplying  (4.40)  by  (4.44)  yields 

RhMh^  = 


(1  +  ^)Ih  +  ahUh 


(4.45a) 


Ultracomputer  Note  143 


Page  31 


where 


0      1 


Uh  = 


Notice  the  similarity  to  (4.4b).    Since 


(4.45b) 


1 
0 


\ltU,\,  =  VV^i"^  ^  ^  . 


a      ,    ah 


we  have 

\RhMh%,l  ^  (1  +  ^  +  -f-)  (4.46) 

which  resembles  (4.8).   Since  U  is  strictly  upper  triangular  it  is  easy  to  compute  MhRh^- 
We  have 

.-1= L 


MhRh 


1  + 


jj^  Uh  -  8A;,)    ,     8  = 


ah 


and 


(Aa)o-  -  1(1  '_  8)0  -'■-  1) 
Since  |l  —  5|  <  1,  we  have  by  direct  computation 


j  >  i    . 


\8A,\, < VfTTT^  , 


which  yields 


lA/./?.-%,i  ^  1  +  Vf  +  1+^    , 


(4.47a) 


(4.47b) 


(4.48) 


(4.49) 


2    '    4    ■     4 
which  resembles  (4.7). 

Example  4.b).    Uniform  discrete  L2  norm  equivalence  of  Mh^  and  Rh^.    Now  consider  the 
inverses.   Multiplying  (4.44)  by  (4.40)  yields 


Mh^Rh  =  (1-  %Ih-  ahUl  +  ah^L  d 


(4.50) 


where  e^  is  as  in  (4.41)  and  v  is  as  in  (4.42).  Notice  the  similarity  to  (4.12).  Here  a  rank 
one  matrix  takes  the  place  of  the  boundary  operator.  The  first  two  terms  in  the  right-hand 
side  are  bounded  as  /i  -  0,  but  the  third  is  not. 
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We  have 


Now 


and 

klo.A  =  h'^   .  (4.51c) 

Thus,  there  exists  K,ho  >  0  such  that  for  h  <  hg 

iMh'Rhkh  ^  Kh-'^  ■  (4.52) 

The  matrix  R^^  M/,  is  messy  to  compute  so  we  only  remark  that  it  can  be  computed  from 
Mh^Rh  by  finding  the  inverse  of  the  first  two  terms  in  (4.50)  and  applying  Lemma  4.3.  The 
result  is  a  bounded  matrix  plus  an  unbounded  rank  one  matrix.   The  net  result  yields 

CoM^h^Rh)  =  \Mh^Rh\o,h  \Rh^Mhkh  ^  Kh-^   ■  (4.53) 

Example  4.c).  Uniform  discrete  H\  norm  equivalence  of  Mh^  and  Rh^.  We  now  consider 
^h^Rh\\,h-  Since  \u\a,h  only  required  the  interior  points  we  could  ignore  the  effect  of  the 
boundary  points.  Here,  however,  we  must  include  (4.37a,c).  If  ii  is  a  grid  function  satisfy- 
ing (4.37a,b,c),  then 

k\i,H  =  J  ("i)2  +  Mlh  +  J  (  2^^)2(w„)2  (4.54) 

where  u.  represents  the  interior  values.  Now  suppose  a  satisfies  (4.37a,b,c)  for  R  (i.e., 
ot  =  <2,7  =  0)  and  i£  satisfies  (4.37a,b,c)  for  M  (i.e.,  a  =  0,  7  =  0).   Then  we  have 

wo  =  0    ,  (4.55a) 

m  =  Mh^RhU.   ,  (4.55b) 

Wn  +  i  =  w„   .  (4.55c) 

Then,  from  (4.54)  we  have 

\iii\lh  =  J  kiP  +  iMh^Rhullh    .  (4.56) 


From  (4.50)  we  have 


Lemma  4.2  yields 


/,        ah  ^       I      , 
wi  =  (1 2~)"i  +  ^hu„    . 
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Uil  :£  ImiI  +  ah\ii\i,h    . 
and 

-I  kiP  ^  -^  ImiP  +  (la  +  a^h)  \u.\Ih  (4.57) 

From  (4.50)  and  Lemma  4.2  we  also  have 

\Mh^Rhll\i,h  ^  \u.\i,h  +  a^U'^a\i,h  +  a\u„\  ]hy.\i,h 

2  lull,/,  +  alilo.A  +  al"nl  (4.58) 

<  (1  +  a(l  +  ^))  \u.\i,h    . 

Combining  (4.57)  and  (4.58)  in  (4.56)  yields 

Imli.A  ^  Ki  \u.\i,h  (4.59) 

for  some  Ki  >  0.    Finally  using  (4.59)  and  Lemma  4.2  we  have 

y.\},h  =  \nL^,h  +  yii:,h 

s  (1  +  ■^)  yb  2(1  +  ^)/^?  kib  (4.60) 

for  some  K2  >  0.    Thus,  Ia/a"^/?aIi,a  is  bounded.    Bounds  on  \Rh^Mh\i,h  can  be  found  in  the 
same  fashion. 

Example  5.    This  example  shows  that  ]NhRh%,h  is  not  uniformly  bounded,  but  ]Nh^Rh\o,h 
and  \Rh^I^h\o,h  are  uniformly  bounded.    It  also  shows  that  ]Nh^Rh\i,h  is  uniformly  bounded. 
Example  5. a).    Uniform  discrete  L2  norm  equivalence  of  N^  and  Rh.    Multiplying  (4.42b)  by 
Rh  we  have 

RhNh^  =  RhMh^  -  yRn^^   .     7  =  ^^ 

Using  (4.44)  we  see 


1  +  a-f-     •  (-^-^l) 


/?A£  =  ahi+  (1  -  -^)e^    ,  (4.62) 

where  1  =  (1,1 1)^.    Combining  (4.45)  and  (4.62)  with  (4.61)  yields 

RhNh^  =  (1  +  -f-)//,  +  ahUn ^^       ,    i  n^ ^  o^y"      (4.63) 

l  +  a  +  -^  \  +  a  -  ^ 

Notice  the  similarity  to  (4.17b).    The  first  three  terms  on  the  right-hand  side  of  (4.63) 
correspond  to   the  terms  in   (4.17b).     We  have  previously  shown   that  the  first  two  are 
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bounded. 
Consider 

.4 


h^liL^}^  ^-n^^l   .  (4.64) 

The  final  term  is  unbounded.   Consider 

^!!5^  -  k,h,  kh.  ~  ^  *-»  .  (4.65) 

Thus,  RN~^  is  bounded  wherever  it  is  defined,  but  the  discrete  analogue  Rf,Ni^^  is  not  uni- 
formly bounded  on  R". 
Example  5.b).    Uniform  discrete  L2  norm  equivalence  ofN^^  and  Rh^. 

Multiplying  (4.42b)  on  the  right  by  Rh  yields 

ah 

1  +  a-f-     '  (4.66) 


Nh^Rh  =  Mh^Rh  -  7£  ii^Rh    ,     y  = 


where 


aJ'Rh  =    -  ahl^  +  (1  +  a  -  -y-)  ej  (4.67) 

Combining  (4.67)  and  (4.50)  into  (4.66)  we  see  the  final  terms  miraculously  cancel  to  yield. 

A.,-^,  =  (1  -  f )  Z.  -  ahUl  +  -f^^  ^'  (4.68) 

I  +  a  —  -~- 

Again,  notice  the  similarity  to  (4.20).  To  show  that  the  last  term  is  bounded  we  have  as  in 
(4.64) 

\h'^^n}=  \h^ly^\}^\    .  (4.69) 

Thus,  combining  (4.46)  and  (4.69)  yields 

K-/?.U  -  (1  +  ^  +  f )  +  VI  (--^^F)  (4.70) 

1  +a  -  — 

The  uniform  boundedness  of  \Rh'^Nh\o,h  can  be  obtained  from  Nh^Rh  by  inverting  the  first 
two  terms  of  (4.68)  and  applying  lemma  4.3.  The  result  is  uniformly  a  bounded  operator 
much  like  (4.47a,b)  and  a  uniformly  bounded  rank  one  operator. 
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Example  5.c).  Uniform  discrete  Hi  norm  equivalence  ofNh^  and  Rh^.  As  in  Example  4.c  the 
first  two  terms  of  (4.68)  are  bounded  in  discrete  Hi  norm.  Since  l^iili,^  =  1,  we  have  by 
lemma  4.2 

]h^iLL^u.\i,h  ^  Ml,;,  i  \uj\h  ^  ]hA,h  k\o,h  ^  ^  \u.\i,h   ■  (4.71) 

Thus,  using  (4.69),  (4.71),  and  Lemma  4.2  we  have 

^   3*  l"B,A  +  Y  \u\\,h  (4.72) 

^  \u\Ih    . 

The  boundedness  of  ]Rh^l^h\\,h  could  be  obtained  in  a  similar  manner  as  Rh^Nh-    We  omit 

this  lengthy  calculation. 

Example  6.    The  matrices  Lh,  Mh,  and  A^;,  are  all  symmetric  positive  definite  with  domains 

D(L)  #  D(M)  #  D(N)  ¥"  D(L).     Thus,   no   two   are   uniformly   equivalent   in    discrete   L2 

norm,  nor  are  their  inverses.    However,  we  have  FoiL)  #  TqIM)  =  ro(A^).    This  example 

will  show  that  Mh  and  Nh  are  uniformly  equivalent  in  discrete  Hi  norm  while  L;,  and  Mh  are 

not.    Finally,  we  will  show  that  Mh  and  Nh  are  uniformly,  spectrally  equivalent  while  Lh  and 

Mh  are  not. 

Example  6.a).     Uniform  discrete  Hi  norm  equivalence  of  N^^  and  Mh^.    Equation  (4.41b) 

yields 

ah 


Mh^Nh  =  I  +  p££l   .     P  = 


1  + 


2 
and  (4.42b)  yields 


ah      '  (4.73a) 


Nh^Mh  =  I  -  yaeJi   ,    7  = 


ah 


(l  +  a-f)     ■  (^-^^b) 


Equations  (4.73a,b)  mimic  (4.23)  and  (4.22)  respectively.    We  have 

h^-Ur  I       ^    \Mh^Nhe^\o,n    ^  |„|    ll-lo,A         ,     ~      a      ,  -u,  ,a-,a^ 

^'  ^'^'^  -         U,H         -  'P'  li:n^  -  ^    =  V?  ^    '  (4.74a) 

Since  MhNjT^  =  (Nh^Mh)*    ,   NhMh^  =  (Mh^Nh)*  we  also  have 

\MhNh%,h  =  V^^f^-^^"''   >  (4.75a) 
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To  find  bounds  in  the  discrete  Hi  norm  notice  that 

l/ii:eIiili,A  =  Wnli.A  l"nl  ^  luli.A 

and  we  see  both  Ia/^WIi_a  and  V^h^Mh\\,h  are  bounded. 

Example  6.b).    Uniform  discrete  Z-2  norm  and  Hi  norm  equivalence  of  Mh^  and  Lh^. 

Equation  (4.41a)  and  (4.42a)  yield 

Mh^Lh  =  I  +  iLel  , 

Lh^Mh  =  I  -  ha  el  . 
As  in  the  previous  examples,  we  see  that 


(4.75b) 


(4.76) 


-3/2 


\Lh'M,]o,h  =  WhLh-%,h  ^  ^  -^^VJ^"" 


(4.77a) 
(4.77b) 

(4.78a) 
(4.78b) 


Equation  (4.76)  implies  that  ]Lh  ^Mh\i,h  is  bounded.    However,  consider  (4.77a).    The  rank 
one  term  satisfies 


Thus,  there  exists  K,ho  >  0  such  that  for  h  <  ho 


(4.79) 


Example  6.c).     Uniform  Spectral  equivalence  of  L,  M  and  N .    Consider  Eqs.  (4.73a)  and 
(4.77a).   The  eigenvalues  of  Mh^Nh  are  1  and 


(1  +  ^ely.)  = 


1  + 


n 


,    ,     ah_    n  +  I 

^  ^    2 


while  the  eigenvalues  of  M/,  ^L^  are  1  and  (1  +  ej£)  =  {I  +  n)  =  h    ^   Thus, 

<NhU.,lL> 


1  < 


<MhU.,lL> 


1  +  a 


while 


1<     <^h^^^>     <;,-! 
<MhlL,lL> 


(4.79a) 


(4.79b) 


and  the  upper  bound  is  achieved. 
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Special  notice  should  be  taken  of  Mh  ^Nh-  Its  spectrum  is  uniformly  bounded  but  its 
discrete  L2  norm  is  not. 

5.   Discretizations 

In  this  section  we  discuss  the  implications  of  the  results  of  Section  3  for  the  families  of 
finite-dimensional  operators  {Ah},  {Bh}  which  arise  from  the  discretization  of  the  invertible, 
uniformly  elliptic  operators  A  and  B .  While  most  of  our  results  apply  to  general  discretiza- 
tion processes,  it  is  convenient  to  present  our  discussion  within  the  finite-element  frame- 
work. 

We  suppose  we  have  a  family  of  finite-dimensional  function  spaces  [Shi  indexed  by  a 
discretization  parameter  /z  -  0  and  invertible  operators  {Ah}  and  {Bh}.  That  is,  for  every 
/  ^  L2(n)  we  have  w,v  ^  H2i^)  the  solutions  of 

Au=f,    Bv  =f  (5.1a) 

and  Uh  and  v^  the  solutions  of  the 

AhUh  =  IPhf   ,    BhUh  =  Phf  (5.1b) 

where  P^  is  a  projection  of  Lii^)  onto  Sh-   With  a  slight  abuse  of  notation  we  write 

Uh  =  Ah^f   ,     Vh  =  BhH   ■  (5.1c) 

We  begin  with  the  theorem  of  Bramble  and  Pasciak  [BP]  which  deals  with  the  uniform 
L2  norm  equivalence  of  {A^'}  and  (5^"^},  given  that  A~^  and  B~^  are  Li  norm  equivalent. 
Then,  we  deal  with  the  corresponding  result  for  the  uniform  L2  norm  equivalence  of  {Ah} 
and  {Bh}.  For  both  of  these  results  we  assume  that  the  discrete  operators  yield  optimal  order 
convergence  estimates  and  satisfy  a  uniform  inverse  condition. 

Definition:  the  family  {Ah}  of  discrete  operators  approximating  the  invertible  uniformly 
elliptic  operator  A  satisfies  Condition  OP  if  there  exists  a  constant  Mi(A),  depending  on  A 
but  not  h,  such  that;  for  every  /  i  Lzi^)  we  have 

\Ah^f  -  A-^f\i^  <  h^MM)\fK   ■  (5.2) 

Definition:  The  family  {Ah}  satisfies  a  uniform  inverse  condition.  Condition  INV  if  there 
exists  a  constant  MiiA),  depending  only  on  A  and  not  on  h,  such  that;  for  every  Vh  ^  Sh  we 
have 

\AhVh\L,^  M2{A)h-'^\vh\L,   .  (5.3) 

Theorem  5.1.  (Bramble  and  Pasciak):  Let  A  and  B  be  two  invertible,  uniformly  elliptic 
operators  and  let  A"^  and  B~^  be  Lj  norm  equivalent.  That  is,  (1.14)  holds.  Specifically, 
we  assume  that  (3.5a,b)  holds.  Let  {Ah}  and  {Bh}  be  the  associated  discretizations.  Let  {Ah} 
and  {5^}  both  satisfy  condition  OP  and  condition  INV.    Then,  there  is  a  constant  M-i{A:B) 
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such  that;  if  {uh}  and  {v„}.    satisfy 

AhUh  =  BhVh  (5.4a) 

then  they  also  satisfy 

WIl,  ^  Mi{A:B)  \vh\l,   ,  (5.4b) 

W\l,  ^  M3(A:B)  \uh\L,  .  (5.4c) 

Proof.  Since  this  result  is  stated  but  not  proven  in  [BP];  we  present  an  elementary  proof. 
Let 

AhUh  =  BhVh  =  fh  (5.5a) 

Uh  =  A-%    ,     n  =  S-V,  (5.5b) 

Then  from  (5.2)  we  have 

\uh\L,  ^  \uh\l,  +  h^MiiA)  VhIl,  . 

From  (3.5a)  and  (5.3)  we  have 

W\l,  ^  Ki(A')K2(B*)  {vhl,  +  h^Mi(A)M2(B)h-^  kL,   .  (5.6) 

However,  we  also  have  from  (5.2) 

\h\L,  ^  kt  +  h^M,(B)  Ifnl,  , 

which  becomes 

kL,  ^  kL,  +  h^MiiB)M2iB)h-^  kl,   .  (5.7) 

substituting  (5.7)  into  (5.6)  we  obtain 

\uh\L,  ^  [Ki(A*)K2(B*)il  +  Mi(B)M2(B))  +  A/i(A)A/2(S)]  kL,   •  (5.8) 

Thus  we  have  obtained  (5.4b).  Interchanging  the  roles  of  u  and  v  yields  (5.4c)  and  the 
theorem  is  proven.    ■ 

Theorem  5.2:  Let  A  and  B  be  invertible,  uniformly  elliptic  operators  which  are  Z,2  norm 
equivalent,  that  is,  (3.2a)  and  (3.2b)  hold.  Assume  that  the  approximating  discrete  opera- 
tors {Ah}  and  {Bh}  satisfy  condition  OP  and  condition  INV.  There  is  a  constant  M4(A:B)  such 
that,  for  every  fh  €  Sh  we  have 

\AhBh^fh\L,  ^  MM-.B)  \fh\L,   ,  (5.9a) 

\BhAh^fh\L,  ^  M4(A:B)  Ifhli,   •  (5.9b) 
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Thus,  the  families  {A^}  and  {5^}  are  uniformly  Li  norm  equivalent. 
Proof.   Let  \fh\L^  =1-    Set 

BkAh^fH  =  gh   .  (5.10) 

We  will  show 

khV,  ^  K2iB)Ki(A)  [1  +  M2(S)Mi(S)]  +  Mi(A)M2(S)    .  (5.11) 

Let 

BA-%  =  g    ,  (5.12a) 

A-'h  =  <\fh    .  (5.12b) 

Ah'fh  =  4)^    .  (5.12c) 

Then 

A^\lH  =fh,     B^\>H  =  g    ,  (5.13a) 

AAh=fh   ,    Bh^H  =  gh    .  (5.13b) 

Let  Wh  ^  S),  be  the  solution  of 

Bhyvh  =  g  (5.14) 

Then,  (3.2b)  and  condition  OP  yields 

^h  -  ^hV,  ^  h^MiiA)  \fhV,  =  h^Mi(A)  (5.15a) 

U>h  -  wJl,  ^  h^Mi(B)\g\L^  <  h^MiiB)K2(B)Ki{A)    .  (5.15b) 

Using  (5.3)  we  have 

Ij  -  ghk,  =  I5a(wa  -  <i>h)\L,  ^  h-'^M2(B)  Uh  -  <^h\L,   ■ 

Using  (5.15a,b)  we  have 

\g  -  ghli,  ^  Mi{B)M2{B)K2{B)Ki{A)  +  Mi(A)M2{B) 

Hence 

\gh\L,  ^  \g\L,  +  \g-  gh\L,  ^  K2(B)KM)[l  +  MiiB)M2(B)]  +  Mi(A)M2(S) 

Thus,  (5.11)  has  been  proved  and  (5.9b)  is  established.    Reversing  the  roles  of  Ah  and  S/, 
completes  the  proof.     ■ 

Theorem  3.1  states  that  the  operators  A~^  and  B~^  are  L2  norm  equivalent  if  and  only 
if  D(A')  =  D{B*).  If  this  condition  is  not  satisfied  then  A^^  and  Bh^  cannot  be  uniformly 
equivalent  in  L2  norm  (see  Theorem  2.9  of  [FMP]).    Our  next  result  deals  with  the  rate  of 
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growth  of  CiJiAh  ^Bf,)  in  the  case 

DiA*)  St  D(5*)    .  (5.16) 

Theorem  5.3.  Let  A  and  B  be  two  invertible,  uniformly  elliptic  operators  which  satisfy 
(5.16).  That  is,  they  do  nol  satisfy  the  hypothesis  of  Theorem  5.1.  Let  {Ah)  and  {S^}  satisfy 
condition  OP .   We  consider  two  cases. 

Case  1.  There  is  a  segment  Yt,  C  Y\{A)  =  Ti{B)  and  on  this  segment  the  boundary  coeffi- 
cients ao,  Po.oti.  Pi  are  not  equal; 

lao  -  3ol  +  Icti  -  Pil  >  0   on  Fs    .  (5.17) 

Then,  there  is  a  constant  ci  >  0  such  that 

\Bh^AH\L,^  cih-'^  ,  (5.18a) 

\Ah^Bh\L,^  ci/i-^   ,  (5.18b) 

and 

CLjiAh^Bh)  ^c}h-^   .  (5.18c) 

Case  2.  There  is  a  segment  T^  C  Tq^A)  D  ri(S).  Then,  there  are  constants  ci,C2  >  0  such 
that 

{Ah^BhIl,^  cih-'^   ,  (5.19a) 

\Bh^Ah\L,^  C2h-^''^   ,  (5.19b) 

and 

CiJiAH^Bh)  ^  C2Cih-^   .  (5.19c) 

Proof.  Consider  Ah'^Bf,  in  both  Case  1  and  Case  2.  The  construction  of  Lemma  3.3  and  the 
discussion  in  Lemma  3.4  show  that,  after  a  smooth  mapping  of  fi  onto  O',  a  domain  which 
satisfies  Condition  ft,  there  is  a  function  Uh  €  DiB)  and  two  constants  c  >  0,  ci  >  0  such 
that 

luhli,  =  1    ,     \Buh\L,  ^  ch-^   .  (5.20a) 

And,  if  Vh  =  A'^Buh  then 

IvaIlj  =  lA"^fluJi.j  >  ci/z"^   .  (5.20b) 

Let  Uh  €  S;,  be  the  solution  of 

BhUh  =  PhBuh  (5.21) 
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while  Vf,  €  Sh  the  solution  of 

AhVh  =  PhBuh  =  BhUh    .  (5.22a) 

That  is 

VA  =  Ah^BhUh    .  (5.22b) 

Then  using  Condition  OP  (5.2) 

W\l,  ^  \uh\L,  +  \uh  -  UhV,  s  1  +  Mi{B)h^\Buh\L,   ■ 
Using  (5.20a)  we  obtain 

\uh\L,  ^  [1  +  Mi(B)F]    .  (5.23) 

On  the  other  hand 

ci/i-^  <  Ivhl,  ^  kL,  +  \n  -  vaL,  ^  \vh\L,  +  h^Mi(A)\Avh\L,   . 
Since  Av),  =  Buh  we  have 

Fi/z-^  <  kt^  +  MiiA)c    .  (5.24) 

Comparing  (5.23)  and  (5.24)  yields  (5.18b).   Reversing  the  roles  of  A/,  and  Bh  yields  (5.18a) 
and  (5.19a). 

A  similar  argument  based  on  the  function  ^h  constructed  in  Lemma  3.3  and  the  discus- 
sion in  Lemma  3.4  yields  (5.19b).   The  Theorem  now  follows.     ■ 

We  now  examine  the  rate  of  growth  of  Ci^{AhBh^)  in  the  case  D{A)  i^  D(B). 
Theorem  5.4.    Let  A  and  B  be  two  invertible  uniformly  elliptic  operations  which  are  not  L2 
norm  equivalent. 

That  is 

D{A)  #  DiB)  (5.25) 

Let  {Ah}  and  {Bh}  satisfy  Condition  OP  with  respect  to  A*  and  S*.    As  in  Theorem  5.3,  we 
consider  two  cases. 

Case  1.    There  is  a  segment  F^  C  TiiA)  =  FiiB)  and  on  this  segment  the  boundary  coeffi- 
cients ao,  Po.  «!.  Pi  are  not  equal.   That  is 

lao  -  Pol  +  lai  -  Pil  >  0  on  Tj    .  (5.26) 

Then,  there  is  a  constant  ci  >  0  such  that 

UaSa"!,  Sci/i-^   ,  (5.27a) 

\BhAh\^cih-'^   ,  (5.27b) 
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and 

Ci^iAhBh^)  ^clh-'^   .  (5.27c) 

Case  2.   There  is  a  segment  Yt,  C  Tq{A)  D  Vi{B).   Then,  there  are  constants  ci,C2  >  0  such 
that 

\BhAh\^  cih-"^  (5.28a) 

\AhBh\  ^  C2/2-3/2  (5.28b) 

and 

ClJ^AhBh^)  2:  C1C2/2-2    .  (5.28c) 

Proof.    We  observe  that  ro(A)  =  ToiA*)  and  ToiB)  =  ToiB*).    The  theorem  then  follows 
from  Theorem  5.3  and  the  observation 

kBirWhi^  =  UhBkX  ,  \iA;)-^B*h\L,  =  \BhAh\  .  ■ 

These  negative  results  immediately  raise  the  question  of  upper  bounds  on  the  rate  at 
which  CiJiBh^Ah)  -^o.    The  general  case  is  not  resolved.    However,  there  is  an  important 

class  of  problems  for  which  we  can  show  that  (5.18)  and  (5.27)  provide  the  exact  asymptotic 
rate. 

As  we  have  said  earlier,  the  theory  of  numerical  methods  for  elliptic  boundary-value 
problem  for  general  boundary  conditions  (1.1b),  (1.1c)  is  rather  sparse.  On  the  other  hand, 
the  finite-element  theory  for  operators  of  Class  N  rests  on  a  solid  basis  and  is  well  under- 
stood. Unfortunately,  even  when  comparing  two  operators  of  Class  N  whose  conormal 
derivatives  are  essentially  different,  our  method  of  proof  requires  that  we  deal  with  the  gen- 
eral case  (1.1c).  Since  we  are  unwilling  to  make  assumptions  (which  may  be  unrealistic) 
about  optimal  order  convergence  for  such  problems-and  hence  complete  transitivity  of  uni- 
form discrete  L2  norm  equivalence  of  Ah^  and  5^^  we  limit  outselves  to  a  special  class  of 
problems.  Fortunately,  this  class  is  rich  enough  to  include  some  cases  of  great  interest,  e.g., 
preconditioning  by  the  discretizations  Bh  of  an  operator  B  with  the  same  leading  part  as  A, 
and  the  case  where  the  leading  parts  of  both  A  and  B  are  diffusion  operators. 

Let  A  be  an  invertible  uniformly  elliptic  operator  of  Class  N  with  coefficients 
^n,  '312  '^22.  ^1.  ^2.  <2o  and  boundary  conditions 

w  =  0  on  ro(A)    ,     jjp  =  au  on  TiiA)  (5.29) 
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Let  'i(x,y)  be  a  smooth  function  which  satisfies 

0  <  70  ^  yix,y)  ^  71   ,    U.y)  ^  ^  (5.30a) 

with  certain  constants  -yo.  7i-    Let  B  be  an  invertible  uniformly  elliptic  operator  of  Class  N 
with  coefficients  bu,  bu,  ^22.  ^l.  *2.  ^o  which  satisfy 

fell  =  7^11    ,     bi2  =  7<3i2   .     ^22  =  7^22   ■  (5.30b) 

Let  To(A)  =  To{B)  and  let  the  boundary  conditions  for  B  be 

M  =  0  on  ro(5)  =  ro(A)    ,     -^  =  pu  on  Ti{B)  =  FiCA)    .  (5.31) 

Let  AT  >  0  and  define  the  uniformly  elliptic  operators  A^,  5^  as  follow.   Let 
A*^M  =    -  {(aiiM;t);c  +  iai2Ux)y  +  (ai2Uy)x  +  (a22Uy)y}  +  Ku    , 

u  =  0  on  ro(A)    , 

■^^  =  au  on  ri(A)    , 

where  a*  is  the  function  associated  with  the  boundary  conditions  for  A*.    Let 
B^v  =    -  {(anvx)x  +  (ai2V;c)y  +  (ai2Vy)x  +  ia22Vy)y}  +  Kv   , 

u  =  0  on  ro(A)  =  ro(S)    , 

^  =  3*"  on  ri(S)    .  (5.33c) 

where  p'  is  the  function  associated  with  the  boundary  conditions  for  B* .   Observe  that 

du  dii  dit  du 

=  7 


(5.32a) 

(5.32b) 

(3.32c) 

(5.33a) 

(5.33b) 

So  that  A^  and  B^  are  self-adjoint  elliptic  operators  of  Class  N  with 

DiA^^)  =  D(A')    ,     D{B^)  =  D{B*)    . 

Let  {Sh)  be  a  family  of  finite-element  spaces  (say  piecewise  polynomial  spaces).    Let 
the  discrete  operators.  Ah,  B/,  be  determined  by  the  bilinear  forms  d(-,-).  b{- ,•).   That  is, 

A)     Consider    the    problem    Au  =  f,f  ^  ^2(^)-     A    function    Uh  ^  Sh  (~^  Hi(D.,rQ)    is    the 
finite-element  approximant  of  the  solution,  u,  if;  for  every  Wh  ^  Sh  C)  Hi(Cl,ro),  Uh 
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satisfies 

aLiuh,Wh)  +  Jf  Wh[aiiuh)x  +  a2(«/.)y  +  aoUh]dxdy 

(5.34) 

-  r  auhWhda  =  (/,Wh)    . 

where  a£.(-,-)  is  given  by  (2.2c). 

B)      Consider  the  problem  Bv  =  f,f  €  LiiH).    A  function  Vh  ^  Sh  (^  Hi(£l,To)  is  the  finite- 
element  approximant  of  the  solution,  v;  if  for  every  w;,  €  5^  fl  HiiCl,rQ),Vf,  satisfies 

bdvh.Wh)  +  n  Wh[biivh)x  +  b2ivh)y  +  boVh]dxdy 

(5.35) 

-  J     ^Vf,Whd(T    =    (J,Wh)      . 

The  operators  A^  and  B^  are  defined  in  a  similar  fashion. 

A.K) 

Consider  the  problem  A^u^  =  f,  f  ^  ^2(^)-  A  function  u^  ^  5^  fl  Hii^lJo)  is  the 
finite-element  approximant  of  the  solution,  u*^,  if;  for  every  w;,  €  S^  n  HiiD.,To),u^ 
satisfies 

adu^,Wh)  +  K{u^,yvh)  -  J  a'aj^w,,  =  (f,Wh)    • 


I  '^'"^ 


B.K) Consider  the  problem  5*^v^  =  f,  f  ^  ^2(^)-  A  function  v^  €  Sh  r^  HiiUJo),  is  the 
finite-element  approximant  of  the  solution,  v*,  if;  for  every  w^  €  S^  H  f/i(ft,ro),v^ 
satisfies 

1 


advLwh)    +    ^(v£wa)    -    I   J   P'v^W;,)    =    (f,WH) 


Theorem  5.5.  Let  A,  B,  A^ ,  B^ ,  Ah,  Bh,  Af,  Bj^  be  as  described  above  and  assume  the 
hypotheses  of  Theorem  3.1.  Assume  there  is  an  ho  >  0  and  a.  Kq  >  0  such  that  for  all 
0  <  h  ^  ho  and  all  K  ^  Ko  the  operators  A^ ,  B^  are  invertible  and  the  discrete  operators 
{Ah},  {Bh},  {Af},  {Bjf}  satisfy  Condition  OP  and  Condition  INV.  Further,  assume  that  the 
finite-element  spaces  5^  satisfy  a  uniform  inverse  condition:  there  is  a  constant  M  >  0  such 
that,  for  all  Wh  ^  Sh  we  have 

IVwhli,  ^  Mh-%h\L,  ■  (5.36) 
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Suppose  there  is  a  segment  Vj  C  FiiA)  and 

a*  =  i  p*on  ri(A)/r3   ,  (5.37a) 

la*  -  —  p1>  Oon  Ts   .  (5.37b) 

Then  there  are  constants  ci,  C2   >  0  such  that 

cih-'^^\Ah^Bh\L,^C2h-'*  ,    h^ho,  (5.38a) 

ci/!-^  :s  \Bh^Ah\L,  ^  C2h-'^  .    h  ^  ho   .  (5.38b) 

Proof.  The  lower  bounds  of  (5.38a),  (5.38b)  were  established  in  Theorem  5.3.  In  fact, 
using  Theorem  5.1  and  the  fact  that  we  may  interchange  the  roles  of  A  and  B,  it  suffices  to 
prove  a  bound  of  the  form 

for  some  K  >  Kq.   Let  K  >  Kq  and  assume  that 

uh.vH  €  Sh  n  f/i(n.ro) 

satisfy 

AUh  =  B^Vh    ,  (5.39a) 

kkj  =  1    •  (5.39b) 

Let  Wh  =  Uh  —  Vh-   Then 

'  3ft 

(5.40) 

+  la-  -  i-p'U   f  \^h\\uh\da    . 
'  aft 

Let  c  =  2(|a*U  +  ly  p'U).   Then  using  (2.1)  on  the  left  and  (2.17)  on  the  right,  we  have 


HA)  |VwJ2,  +  K  \wh\l^  <  ^0  [IVvvaL^  •  kli,  +  klZj 


(5.41) 


Apply  (3.46)  with 


To  obtain 


€=     AMI  r.^     _£^ 

^0      '  HA) 
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with  ci  =  cl(/min(-^^    .   y  ^:).   Using  (5.36)  we  have 

Thus,  using  the  triangle  inequality  we  have 

Mr^A^unl,  =  kli,  ^  IuJl,  +  \uH  -  v^Il,   . 
That  is 

The  theorem  now  follows.     ■ 

Observe  that  this  proof  remains  valid  if  we  only  require  the  equalities  (5.30)  to  hold  on 
ri(A). 

Remark:  While  this  theorem  seems  to  allow  the  case  ro(A)  #  0  and  TiiA)  #  0,  H2  esti- 
mates are  essential  to  (i)  reduce  the  general  problem  to  a  discussion  of  A^  and  B^ ,  (ii)  apply 
Theorem  5.3  in  order  to  obtain  the  lower  bounds  of  (5.38).  However,  consider  the  special 
case  where  the  leading  part  of  A  =  leading  part  of  S.  In  this  case,  the  argument  of  Lemma 
3.6  enables  us  to  handle  the  case  where  ao  —  ^  or  ^0  —  K.  The  argument  of  Theorem  5.5 
then  proves  the  upper  bounds  of  (5.38). 

We  close  this  section  with  two  results  on  /fi-norm  equivalence  of  {A^^}  and  {8^^}. 
Theorem  5.6.   Let  A  and  B  be  two  invertible  uniformly  elliptic  operators  of  Class  N  with 

To(A)  =  To(B)    . 

Let  {Aft}  and  {S^}  be  the  associated  discretizations.  Suppose  these  discretizations  are  of  the 
form  (5.34),  (5.35);  that  is,  the  discretization  is  obtained  by  merely  restricting  (2.14)  to  the 
subspace  Sh  H  Hi(0,,ro(A)).  Let  there  be  an  /iq  >  0  such  that,  for  all  0  <  ^  ^  ho.  Ah  and 
Bh  are  uniformly  invertible  on  HliilXo),  the  space  of  bounded  linear  functionals,  on 
Hi(n,To).   Then 

\(AH^BH)\H,^nd\Bh^AH\H, 

are  uniformly  bounded  for  0  <  /i  :S  /zq. 

Proof.   The  arguments  in  the  proof  of  Theorem  3.2  go  over  essentially  word-for-word.     ■ 
Theorem  5.7.    Let  A  and  B  be  two  invertible  uniformly  elliptic  operators.    Let  the  discrete 
operators  satisfy  the  Hi  analog  of  condition  OP.   That  is,  there  are  constants  MsiA),  MsiB) 
such  that,  for  all/  €  L2 

Uft-y  -  A -1/1//^  <  hM5(A)\f\L^  (5.42) 

with  a  similar  estimate  for  B~^  and  Bh^.  Suppose  ToiA)  n  TiiB)  ¥=  0.  Then  there  is  a  con- 
stant c  >  0  such  that 
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Iflft-Uj;/,  s  a-^   .  (5.43) 

Proof.    We  will  prove  the  theorem  in  the  case  described  in  Lemma  3.4,  Case  2.    We  recall 
that  proof.   Using  the  function  4>t,  given  by  (3.13b)  we  construct  iI/t,  so  that 

and 

W^   =    Vji^   -    Cj)^    =   T1-3/2h;i     ,       0  <  /I    ^    1     , 

where  vvi  is  a  fixed  function. 

Unfortunately,  V^^^u^  is  of  order  "n"^.    However,  we  may  consider  q^ 

q^  =  -^^r,   .  (5.44) 


Then,  (3.10c)  asserts  that 


Let 


Then 


and 


A  =  ^n  'I'll  • 


fy\  =  q-n  +  Tl'^i    .  (5.45) 
where  wj  is  a  fixed  function.   Let  n  =  h  m  (5.45)  and  let  ma.v^  satisfy 

AhUh  =  Aqh  (5.46a) 

BhVh  =  Aqh  =  Bfh    .  (5.46b) 
Then, 

l"J«,  ^  hhiu,  +  \qh-  Uh\H,  (5.47a) 

\fh\H,  ^  kl//,  +  k  -  fh\H,   ■  (5.47b) 
Using  (5.42),  (1.9),  (5.44)and  (3.10c)  we  see  that 

hh  -  wjf/,  s  hMsiA)  \Aqh\L,  ^  hMs{A)\qh\n,  ^  A/5(A)c  ,                    (5.48a) 

k  -  fh\H,  ^  hMsiB)  lAqnl,  ^  hMsiBMn,  ^  M5(B)c  ,                   (5.48b) 
Therefore,  using  (5.47a)  we  see  that 

Ultracomputer  Note  143  Page  48 


WhIh  —  c    ,   some  constant   , 
while  (5.47b),  (5.48b),  and  (5.45)  imply  that 

and  the  theorem  follows.     ■ 

6.   Summary 

In  this  paper,  we  have  discussed  the  effect  of  boundary  conditions  on  'the  norm 
equivalence  of  uniformly  elliptic  operators  and  their  discretizations.  The  motivation  for 
examining  equivalent  operators  stems  from  the  desire  to  construct  uniformly  equivalent  fam- 
ilies of  discrete  approximations.  It  was  shown  in  [FMP]  that  the  discretizations  cannot  be 
uniformly  equivalent  unless  their  limit  operators  are  equivalent. 

We  have  examined  the  Li  norm  equivalence  of  the  forward  operators,  A  and  B  and  the 
L2  norm  and  H\  norm  equivalence  of  the  inverse  operators  A~^  and  B~^.  The  results  on  L2 
norm  equivalence  depend  upon  the  H2  regularity  bounds  (1.7)  and  (1.9).  Equivalence  in  L2 
norm  in  the  absence  of  Hi  regularity  is  an  open  question.  The  results  on  Hi  norm 
equivalence  do  not  depend  on  Hj  regularity  but  are  restricted  to  operators  of  Class  N.  Both 
Li  and  H\  results  apply  to  a  wide  class  of  operators. 

In  brief,  in  the  presence  of  H2  regularity  for  all  relevant  uniformly  elliptic  operators  we 
have  shown  that 

I.  A~^  and  B~^  are  Z-2  norm  equivalent  on  L2(fi)  if  and  only  if  D{A*-)  =  D(B*). 

II.  A  and  B  are  L2  norm  equivalent  on  Lii^)  if  and  only  if  D(A)  =  D(B). 

These  results  show  that  the  conventional  wisdom  of  choosing  B  to  be  the  leading  part 
of  A  with  the  same  boundary  conditions  as  A  is  only  appropriate  if  post  conditioning  is  being 
used,  that  is,  if  ClJ^AB'^)  is  to  be  bounded.    If  CiJiB'^A)  is  to  be  bounded,  which  is  most 

often  the  case  when  preconditioning  is  used,  then  we  must  choose  the  boundary  conditions 
of  B  so  that  D{A*)  =  D(B*). 

The  result  on  Hi  norm  equivalence  can  be  simply  stated. 

III.  A~^  and  B~^  are  Hi  norm  equivalent  on  Lii^)  if  and  only  if  To(A)  =  FoiB). 

Since  H2  regularity  is  not  required  for  this  result  it  is  appropriate  to  discuss  general 
boundary  conditions  for  operators  of  Class  N.  Also,  this  result  is  not  based  upon  the 
requirement  that  A  and  B  yield  positive  definite  bilinear  forms.  However,  in  this  case  little 
theory  exists  on  how  to  construct  discrete  iterations  that  are  based  on  a  discrete  Hi  norm. 

In  Section  5,  we  examined  families  of  discretizations.  We  established  some  rates  of 
growth  for  the  various  discrete  conditions  when  the  limit  operators  did  not  satisfy  I,  II,  or 
III  above.  We  also  showed  that  if  discretizations  satisfy  Condition  OP  and  Condition  INV 
then 
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IV.  {Ah^  and  {5;,  ^}  are  L2  norm  equivalent  on  1-2(0)  if  and  only  if  A    ^  and  B    ^  are  L2 
norm  equivalent  on  Lji^)- 

V.  {Ah}  and  {fi/,}  are  L2  norm  equivalent  on  L2(fi)  if  and  only  if  A  and  B  are  L2  norm 
equivalent  on  Liiil). 

The  first  result  is  an  adaptation  of  a  result  by  [BP]. 

We  also  show  that  A  and  B  are  in  Class  N  and  if  {A/,}  and  {Bf,}  were  derived  from  res- 
tricting the  weak  forms  of  A  and  B  to  the  subspace  Sh  and  are  uniformly  invertible  then 

VI.  [Ah^  and  (5^"^}  are  f/i  norm  equivalent  on  £-2(0)  if  and  only  if  A~^  and  B~^  are  //i 
norm  equivalent. 

Thus  uniform  discrete  L2  equivalence  can  be  easily  established  if  Condition  OP  and 
Condition  INV  hold.  The  latter  is  always  almost  present.  If  Condition  OP  does  not  hold 
then  our  techniques  are  inadequate.  Some  specific  examples  of  finite  difference  operators 
are  examined  in  [FMP]. 

The  discrete  Hi  norm  result  is  rather  general.  It  applies  to  a  wide  class  of  finite  ele- 
ment discretizations  using  the  weak  form  of  the  operators.  Again  we  remark  that  the  theory 
for  discrete  iterations  based  on  the  discrete  H\  norm  is  limited  for  nonself-adjoint  positive 
definite  operators. 

We  also  remark  that  two  different  discretizations  of  the  same  operator  which  act  on  the 
same  subspace  will  yield  uniformly  Lj  norm  equivalent  families  if  they  both  possess  Condi- 
tion OP  and  Condition  INV. 

Finally,  we  emphasize  that,  while  condition  numbers  play  an  important  role  in  conver- 
gence bounds  for  iterative  methods,  they  do  not  tell  the  whole  story.  Firstly,  a  condition 
may  be  bounded  with  respect  to  h,  but  very  large.  Thus,  some  equivalent  operators  may 
provide  a  poor  preconditioning.  Secondly,  preconditioning  by  a  nonequivalent  operator  may 
produce  a  bounded  operator  plus  an  unbounded  operator  of  low  rank.  For  example,  notice 
that  in  each  of  the  examples  in  Section  4,  the  unbounded  portion  was  a  rank  one  operator. 
In  the  proper  context,  conjugate  gradient  methods  would  converge  much  faster  than  bounds 
using  the  condition  number  might  imply.  In  fact,  in  Example  (6.b)  a  conjugate  gradient 
iteration  would  converge  in  two  steps  despite  the  unbounded  condition  number.  In  more 
than  one  dimension  the  situation  is  more  complicated.  The  rank  of  the  unbounded  portion 
appears  to  increase  in  proportion  to  the  square  root  of  the  number  of  grid  points  on  a  boun- 
dary segment  with  conflicting  boundary  conditions.  An  analysis  of  this  phenomenon 
together  with  other  numerical  issues  will  be  discussed  in  a  future  report. 
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